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COMPETITIVE EROSION IS CONFORMALLY INVARIANT 


SHIRSHENDU GANGULY AND YUVAL PERES 


Abstract. We study a graph-theoretic model of interface dynamics called competitive erosion. 
Each vertex of the graph is occupied by a particle, which can be either red or blue. New red and 
blue particles are emitted alternately from their respective bases and perform random walk. On 
encountering a particle of the opposite color they remove it and occupy its position. We consider 
competitive erosion on discretizations of smooth planar simply connected domains. The main 
result of this article shows that at stationarity, with high probability the blue and the red regions 
are separated by an orthogonal circular arc on the disc and by a suitable hyperbolic geodesic on a 
general ‘smooth’ simply connected domain. This establishes conformal invariance of the model. 
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Figure 1. Competitive erosion on the unit disc, spontaneously forms an interface. 
i and —i form the red source and blue source respectively. In the initial state each 
vertex of a mesh of size ^ is independently colored blue with probability 1/3 and 
red otherwise. After 5000 time steps of competitive erosion, the red and blue regions 
have separated with the interface being an orthogonal circular arc. 


1. Interface dynamics 

In 2003, Propp introduced competitive erosion: a graph-theoretic model of interface dynamics 
maintained in equilibrium by equal and opposing forces. The model has the following underlying 
data: 

• A finite connected graph G = (V, E ) with vertex set V and edge set E. 

• Probability measures and /X 2 on V. 

• An integer 0 < k < jfV — 1. 
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Competitive erosion is a discrete-time Markov chain (S(t))t> o on the space 


{5ch: \S\ = k}. 


One time step is defined by 

s(t+i) = (s(t)u{x t })-m} (i.i) 

where X-f is the hrst site in S(t) c visited by a simple random walk whose starting point has distribu¬ 
tion fii ; and 1) is the first site in S(t) U {Xt} visited by an independent simple random walk whose 
starting point has distribution H 2 - As in Fig. 1 we will think of every vertex colored either ‘blue’ or 
‘red’ and S(t) will denote the set of blue vertices. If the distributions p,\ and /r 2 have well-separated 
supports, one expects that these dynamics separate the graph into coherent red and blue territories. 


Thus competitive erosion is a competing version of Internal DLA which is a fundamental growth 
model; first proposed by Meakin and Deutch in 1986 as a model of industrial chemical processes 
such as electropolishing, corrosion and etching. Competing particle systems modeling co-exitence 
of various species etc, have been the subject of intense study in physical sciences as well as mathe¬ 
matics: competing versions of the Richardson model were considered in [2, 18], processes modeling 
annihilation between competing species were studied in a series of works, see for e.g. [ , 8]. The 
study of competitive erosion was initiated in [15] where the reader can hnd a detailed introduction 
of the process. The underlying graph considered in [ ] was the cylinder (the product of a path and 

a cycle). However the original question asked by Propp was in the setting where the underlying 
graph is a discrete approximation to a general smooth simply connected domain. The process was 
predicted to exhibit conformal invariance. Confirming this is the goal of the article. More on 
comparison with Internal DLA is presented in Section 1.4. 


Many 2D lattice models of physical phenomena are believed to have conformally invariant scaling 
limits: percolation, Ising model, self-avoiding polymers, dimer models etc. The last decade and a 
half has witnessed tremendous progress in understanding such models. Below we mention some 
of the seminal works; conformal invariance of various observables related to dimer models, was 
established by Kenyon (see [21, 22, 23]), convergence to conformally invariant scaling limit for loop 
erased random walk was established by Lawler, Schramm, Werner [ !5], conformal invariance in 
percolation and spin systems was proven by Smirnov, see [29] and the references therein. 

1.1. Conformal Invariance and informal setup. We consider the case when the underlying 
graph is a discretization of a smooth planar simply connected domain HJ i.e. V = U n := ^Z 2 n U 
and k = [a:|U n |J for some fixed a G (0,1/2], Consider the case when the underlying domain is 
the unit disc D thought of as a subset of the complex plane C. Simulations (Fig. 1) show if the 
measures p,\ and fi 2 are Dirac measures on —i, i, then after running the process for some time the 
blue and red regions are separated with the boundary being an orthogonal circular arc such that 
the blue region has a fraction of the total area. Let us call this region B/ a \. For a general domain 
HJ with two points (aq, X 2 ) on the boundary, one can obtain the corresponding U( a ) by conformally 
mapping D to U with —i, i mapped to x 1 , aq respectively and using the area measure on U to ensure 
that U( a ) contains a fraction of area(U), (see Fig 3). Precise definitions are provided in subsection 
1.3. 

It was predicted by Propp in 2003 based on the conformally invariant nature of harmonic measure 
that competitive erosion should exhibit conformal invariance. That is, for any reasonably regular 
domain HJ and the erosion chain on U n with k = [a|U n |J and the blue and red sources being x\. aq 
respectively the blue region should “converge” to U( a ) as n grows large. Our main result confirms 
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this. 

However for technical reasons we need Hi,H 2 to be “smooth” and supported in the interior on U 
instead of being point masses. In this article we take them to be uniform measures on all lattice 
points inside small ‘blobs’ of radii 5 /4 lying inside U, and also at distance 5 from the points x\, X 2 
respectively. See Fig 2 i. 

The main result then involves sending 5 to 0 to “approximate” the Dirac measures at X\,X 2 - We 
now present the statement of our main result. 

To improve readability at this point we postpone the precise description of the technical assumptions 
including the definitions of the blobs. These are discussed right after the statement of the theorem 
(subsection 1.3). 



Figure 2. On the left: uniform measure on the two blobs close to the points x\,X 2 
act as “smooth” approximations of the corresponding Dirac measures. On the right: 

The red and blue territories have separated out after running the erosion chain with 
the blue region being “close” to HJ a if the blobs are small enough. 

1.2. Main result. Recall that a € (0, is the fraction of blue vertices. Let a simply connected 
domain “U C C is smooth” mean that the boundary of U is an analytic curve (equivalently the 
conformal map from U to D has a conformal extension across the boundary, see [ , Prop 3.1]). 

Let 

Qe,n := {5CU„: |5| = |a|U„|J, U (Q _ e) nO„c5c U (a+e) n U„}, (1.2) 

Note that because of the smoothness of the domain U, the boundaries of U( Q ) and U( a ± e ) are 0(e) 
away from each other, a will be fixed throughout the rest of the article. 

Theorem 1.1. (Main Result) Let HJ C C be a smooth simply connected domain. Consider the 
competitive erosion chain on U n = U n ^Z 2 with blob radii 5 /4 and Lo:|U n |J blue vertices. Then, 
given e > 0 for 6 < 60 (e) there exists a positive constant d = d(e, 6 ) such that for all large enough 
n = 2 m , 

^,n(Ge,n)>l-e~ dnl/3 (1.3) 

where S is the set of blue vertices sampled from Trs, n , the stationary measure of the chain. 

See Fig 1 and 2 ii. Note that for brevity we suppress the centers of the blobs in the notation 1 : 5 ^. 
The centers are points chosen at distance at least f from the boundary of U and also at distance 
6 , from X\,X 2 respectively (see Fig 4), formally described in Section 1.3. 

Informally the above theorem says the following: consider the erosion chain on the discretiza¬ 
tion of any smooth domain U with points xi,X 2 on the boundary. Suppose that a fraction of the 
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vertices are blue. If the blobs from which the blue and red random walks start are small enough 
and close enough to x\ and then as the mesh size goes to zero, at stationarity, the blue region 
looks like the set HJ( a ), in the sense that on the ‘U( Q ) side’ of a small band around the boundary 
of U (a) , the vertices are all blue and similarly on the ‘other side’, all vertices are red. As discussed 
above, the sets HJ( Q ) for various domains U can be obtained from each other via conformal maps. 
Thus Theorem 1.1 establishes that competitive erosion is conformally invariant. 

Note that the mesh size goes to 0 at dyadic scales. This is for technical convenience. Reasons are 
elaborated in Section 1.5 iii. 


1.3. Formal definitions and setup. In this section we make precise all the definitions and the 
entire setup leading to the statement of Theorem 1.1. We will also present a more quantitative 
version of the result. 

For any two points x, y E C, d(x , y ) will be used to denote the euclidean distance between them. 
For any set A C C and any i£C denote by d(x,A ), the distance between the point and the set. 
D will be used to denote the unit disc centered at the origin in the complex plane. B(x, e) denotes 
the open euclidean ball of radius e with center x. 


Throughout the rest of the article all our domains will be bounded, simply connected, smooth (as 
described before the statement of Theorem 1.1) and hence we will drop the adjectives. 

Definition 1. For any domain B C C denote by dB the boundary of B. 


For a domain HJ and points x\, X2 E dU let 

0:D->U (1.4) 

if : HJ —> D 

be conformal maps such that (fro if, ifr o <fr are the identity maps on the respective domains and 
ifr(xi) = —i, tfr(x2) = i- The existence of such maps is guaranteed by the Riemann Mapping 
Theorem. See for eg: [1, Chapter 6]. In fact there exists more than one pair of ((fr,ifr) since a 
conformal map between domains has three degrees of freedom and here we have fixed the value at 
only two points. However we choose a particular pair (cfr , ifr) assumed to be fixed throughout the 
rest of the article. For any /9sM define 


Ug := jz E U : — log 


ifr(z) — i 
ifr{z) + i 



(1.5) 


The constant — is not important. It falls out of some natural integrals involving the Brownian 
motion heat kernel appearing later in the article. Thus for the disc, is a region containing —i and 
enclosed by an orthogonal circular arc (geodesic with respect to the hyperbolic metric) symmetric 
with respect to i and —i. For a general domain U, we transfer the regions via conformal maps and 
the boundaries still remain geodesics as they are conformally invariant. See Fig. 3. However for 
our purposes we need an area parametrization of the regions Ug. Given a E (0, 1), let /3 = /3(a) be 
such that 


Let 


area(U^) = a area(U). 


U (a) • Ug( Q p 

We now describe the set up in the statement of Theorem 1.1 formally. 


( 1 . 6 ) 

(1.7) 
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Figure 3. Hyperbolic geodesics are circular arcs on the disc. They are invariant 
under conformal maps. However since conformal maps are not area preserving D a 
can get mapped to U Q / for some a' ^ a . 

Setup 1. Given U we take U„ = Un 1 i 2 ), as our vertex set. As the edges of our graph we take 
the usual nearest-neighbor edges of U n though of as a subset of ^Z 2 . However we delete every 
such edge which intersects U c . Since U is assumed to be smooth, U n will be a connected for large 
enough n. See Remark 1.1 below. 

Fix a £ (0, ^], and take k = [o;ITJ^|J. Let X\,X 2 G <9U. The following defines the blobs men¬ 
tioned in the statement of Theorem 1.1. For small enough <5 > 0 let yi, y 2 £ U be such that 

d(xi,yi ) = 6 

d(yi,d U) > 5/2. 

For i = 1, 2, let U ? ; = B(yi, |). As discrete approximations of U* we take 

= i, —) n u n 

where Zi tTl G ^Z 2 is the closest lattice point to y*. Define fii to be the uniform measure on Uj jn . 

Thus the blobs are the sets Ui and U 2 and n /s are uniform over U i >n . Note that in the above, y^’s 
were just required to satisfy certain properties and other than that were completely arbitrary. Also 
we abuse notation a little in the definition of the blobs, n should be thought of as large and hence 
U n (the underlying graph) should not be confused with the blobs Ui and U 2 . 

Remark 1.1. The smoothness assumption on U allows us to choose the y/s. This is formally 
proved in Corollary 9.1. See Fig. 4. The connectedness of U n for large enough n follows since near 
the boundary U looks like a half plane locally, see (9.1). 

It will be instructive to prove first the following weaker version of Theorem 1.1. 

1.3.1. Statistical version of Theorem 1.1. As mentioned already, to highlight the symmetrical roles 
played by S(t) and its complement (see (1.1)), we will think of the states as 2-colorings (blue and 
red) rather than sets. In this language S(t ) will denote the set of all blue vertices. 

Formally define, 




n := Q n := {a G {1, 2} u ", #{x G U n : a{x) = 1} = La|U n |J} 
n' := n' n := {cr G {1, 2} u ", #{x G U n : a{x) = 1} = La|U n |J + 1}. 


( 1 . 8 ) 
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Figure 4. y^s are points at distance <5 from Xj’s. They are at distance at least | 
from (9U. The blobs are discs of radius | centered at the y^s. 


We let, 


a t (x) 


1 x € S(t) 

2 X S(t). 


(1.9) 


Thus the competitive erosion chain S(t), can be thought of as a Markov chain at , t = 0,1,2,..., on 
Cl. However for later purposes we need the following notation which keeps track of the full process, 


at, t = 0,1/2,1,3/2, 2,... (1.10) 

where at 6 14 or at 6 JT depending on whether t is an integer or a half integer. Clearly a single 
time step of the chain consists of a step from Cl to Cl' followed by a step from Cl' back to Cl. 

Given e > 0, define A t := A n)€ to be the set of all configurations a 6 Cl n such that 

#{x G U( Q ) : a(x) =2} < en 2 (1.11) 

#{z € U( q) : a(x) = 1} < en 2 , 

where U( Q ) was dehned in (1.7). Thus informally if e is small, Cl n ^ is the set of all configurations 
where the amount of “dust” particles of the wrong color has small density. See Fig. 5. The next 



Figure 5. A typical configuration in Cl ne for small e with small number of particles 
of the “wrong” color. By Theorem 1.2 the stationary measure of the erosion chain 
concentrates on such configurations. 

result shows that the stationary measure of A n)t is large. 
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Theorem 1.2. . Let U be as in Setup 1 and e > 0. Then there exists do = do(e) such that for all 
5 < S 0 and n = 2 m > N{5), 

7D,n(Ai,e) > 1 “ 

where D = D(e, <5, U) > 0, where ns t n, the stationary measure of the erosion chain on U n with 5 as 
in Setup 1. 

1.4. Comparison with IDLA. Internal diffusion limited aggregation (IDLA) is a fundamental 
model of a random interface moving in a monotone (outward) fashion. IDLA involves only one 
species with an ever-growing territory 

I(t + 1) = I(t) U { X t } 

where Xt is the first site in I(t) c visited by a simple random walk whose starting point has distri¬ 
bution // 1 . Competitive erosion can be viewed as a symmetrized version of IDLA: whereas I(t) and 
I(t) c play asymmetric roles, S(t) and S(t) c play symmetric roles in (1.1). IDLA on a finite graph is 
only defined up to the finite time t when I(t) is the entire vertex set. For this reason, the IDLA is 
usually studied on an infinite graph and the theorems about IDLA are limit theorems: asymptotic 
shape [24], order of the fluctuations [3, 4, 19], and distributional limit of the fluctuations [20]. In 
contrast, competitive erosion on a finite graph is defined for all times, so it is natural to ask about 
its stationary distribution. To appreciate the difference in character between IDLA and competitive 
erosion, note that the stationary distribution of the latter assigns tiny but positive probability to 
configurations that look very different to the final figure in Fig. 1. Thus competitive erosion will 
occasionally form these exceptional configurations. 

1.5. Remarks about Theorems 1.1, 1.2 and Setup 1. 

i. Since U is smooth, using Schwarz reflection cf and hence if can be extended conformally 
across the boundary onto some neighborhoods of D and U. In particular this implies that 
\<ff\ and \if'\ are bounded away from 0 and ooonlJ and D respectively. See [28, Prop 3.1]. 
This bi-Lipschitz nature of the maps will be used in several distortion estimates throughout 
the rest of the article. 

ii. Note that the blob sources U* lie entirely in the interior of HJ and the measures m should 

be thought of as “smooth” approximations to point masses at points x\, X 2 as 5 —> 0. Also 
by choice = |U 2 ,n|- This will be technically convenient. 

iii. Lastly we discuss the choice of the dyadic mesh sizes in the statements of Theorems 1.1 

and 1.2. The technical core of the proof of the theorems involve convergence of the random 
walk on U n under proper scaling of time to Reflected Brownian motion on HJ. The proof of 
this appears in [9] and subsequent local CLT estimates have been obtained in [ ]. However 

both the above papers assume dyadic discretization in their proofs. Since the proof of 
our results rely heavily on the aforementioned convergence results we stick to the dyadic 
discretization in our statements as well. We end by mentioning that it has been pointed 
out (via personal communication) by the authors of [ ] and [ ] that the convergence holds 

true even when the mesh size at the n th step is ^ instead of ^ and the latter was chosen 
for technical convenience. 

2. Sketch of the proofs and organization of the article 

To prove Theorem 1.2 i.e. the stationary distribution concentrates on the small set A e := A n e , we 
identify a Lyapunov function w(-) on the state space. That is a function which attains its global 
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maximum in A e and increases in expectation in one step of the process when starting outside A e , 
i.e. if a q £ A e , 

E(w(<Ji) — w(cto) I o"o) > a > 0. (2-1) 

For more on Lyapunov functions see [ ]. To construct such a function w(-) we proceed by defining 

the following discrete Green function: for any x E U n , 

2 n 2 r°° 

G n (x ) = / P x (X(t) E Ui, n ) - P x (X{t) E U 2 ,„) dt, (2.2) 

I P l,n | Jo 

where P x (-) is the measure induced by the continuous time random walk on the graph U n started 
from x with mean holding time . This choice of holding time is made to ensure that the random 
walk converges to Reflecting Brownian motion on U. Thus the above function is the expected 
difference in the amount of time random walk spends in Ui n and in P^n respectively. Now for any 
(TGllUSl'we define the weight function, 

w{a) ■= ^2 G n (x), 

xGBi(rr) 

where for a E 12 U fL and * = 1,2 

Bi(a) := {x E U n : cr(x) = i}. (2.3) 

Recalling S(t) from the dehnition of competitive erosion we note that S(t) = Now recall 

the map if from (1.4). 

Theorem 5.1 states that: up to translation by a constant G n (x) is close to the function ^ log | 2(x)+l I 
if the blob sizes are small and n is large. 

This justifies why w(cr) is maximized if all the blue vertices are in U( Q ) (see (1.7)). The proof of 
the above is technical and involves developing a convergence theory for discrete Green function 
with Neumann boundary conditions (this is done in [16]). The proofs rely on local central limit 
theorems for convergence of random walk on U n to Reflected Brownian motion in U. Some of these 
estimates were obtained very recently in [9], [11] and [13]. 

The technical core in the paper consists of the following: 

• Proof of Theorem 4.1, which uses an electrical resistance argument involving Rayleigh’s 
monotonicity principle to establish the drift (2.1). 

Once (2.1) is established we use Azuma’s inequality to argue that the process w(at) spends nearly 
all its time in a neighborhood of its maximum. These ingredients together with a general estimate 
relating hitting times to stationary distributions (Lemma 4.7) establish Theorem 1.2. 

To prove Theorem 1.1 from Theorem 1.2 we show that starting erosion from a configuration a E A e , 
the remaining dust particles (red particles in U( a _ e )) get wiped off quickly. This is shown using 
robust IDLA estimates on U n starting from the random environment a. Proof of such estimates 
constitutes the remaining technical challenge. 

2.1. The level set heuristic. We now mention a heuristic that justifies the construction of the 
function w(-) outlined in the previous subsection and predicts the location of the competitive erosion 
interface. Consider the case of “well-separated” measures g.i, g/ 2 - on general finite connected graph, 
which we assume for simplicity to be r-regular. Let g be a function on the vertices satisfying 

Ag = hi - M2 (2.4) 

where A denotes the Laplacian A g{x) := ^ — g{y)) and the sum is over vertices y 

neighboring x. Since the graph is assumed connected, the kernel of A is one-dimensional consisting 
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of the constant functions, so that equation (2.4) determines g up to an additive constant. 

The boundary of a set of vertices S C V is the set 

dS = {v£ V\S : v ~ w for some w G S}, (2.5) 

where v ~ w denotes that they are neighboring vertices. Consider a partition of the vertex set 
V = Si U B U S 2 where dS\ = B = dS- 2 - (Think of Si and S 2 as the blue and red territories 
respectively, of the sort we might expect to see in equilibrium, and the sites in their common 
boundary B have indeterminate color.) 

Let gt for i = 1, 2 be the Green function for random walk started according to gi and stopped on 
exiting Si. These functions satisfy 

A gi = m on Si, 
gi = 0 on Sf. 

The probability that simple random walk started according to Hi first exits S* at x £ B is —A gi(x). 
To maintain equilibrium in competitive erosion, we seek a partition such that Agi ~ Ag 2 he. they 
are roughly equal, on B, that is 

A (01 - c/ 2 ) ~ Mi - M 2 , 

(exact equality holds except on B). Thus by (2.4), the function g — (gi — < 72 ) is approximately 
constant. Since gi vanishes on B, the equilibrium interface B should have the property that 

g is approximately constant on B. 

The partition that comes closest to achieving this goal takes S\ to be the level set 

5*1 = {x : g(x) < K} (2.6) 

for a cutoff K chosen to make #Si = k. An application of the maximum principle shows that for 
this choice of Si, the maximum and minimum values of g — (gi — c/ 2 ) differ by at most 

max \g(x)-g(y)\, 
xeSi,y£Si,x~y 

suggesting that the right notion of “well-separated” measures gi an d M 2 is that the resulting function 
g has small gradient. 

In our case as in Setup 1, gi is the uniform measure on Uj iTl . It is easy to check and shown in 
(5.4) that A G n = Mi — M 2 where G n was defined in (2.2). This shows why the choice of the weight 
function is “natural”. 

2.2. Organization of the paper. In Section 3 we show that the competitive erosion chain has a 
unique communicating class of recurrent states. The proof of the main result appears in Section 7. 
Theorem 1.2 is proven in Section 4. The key idea behind the proof of Theorem 1.2 is the construction 
of a suitable Lyapunov function as discussed in the previous subsection. This construction involves 
defining a suitable potential function. This is done is Section 5. Theorem 4.1 states that the 
constructed function indeed can be used as a Lyapunov function by establishing the drift condition 
(2.1). The proof of this theorem poses the main technical challenge which we address in Section 
6. Theorem 5.1 is one of the main convergence results that we use in this paper which states that 
the Lyapunov function is asymptotically conformally invariant. We discuss a sketch of the proof in 
Section 8. The complete proof appears in [16]. 

Theorem 1.2 follows from Theorem 4.1 and hitting time estimates for sub-martingales and a general 
result in markov chains relation hitting times to stationary measure. Section 4 is devoted to this. 
The proof of the main result relies on additional robust IDLA estimates on HJ n which are provided 
in the Appendix. Some basic geometric facts are proved in Section 9. 
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2.3. Assumptions, notations and conventions. 

i. Through out the article, by random walk on U„, we will mean the continuous time random 
walk with exp(2n 2 ) waiting times unless specifically mentioned otherwise. This is done to 
ensure that the random walk density converges to that of RBM. A fact which would be 
used heavily through out the article. 

We now recall some of the notations already used and introduced some new notations to be used 
throughout the rest of the article. We will denote the complex plane by C. As already used earlier 
for any two points x, y £ C, d(x , y) will denote the euclidean distance between them. Also for any 
set A C C and any x E C denote by d(x,A ), the distance between the point and the set and simi¬ 
larly d(A,B) denotes the distance between two sets A, B C C. B(x,e) denotes the open euclidean 
ball of radius e with center x. 

For any process and a subset A of the corresponding state space r(A) will denote the hitting time 
of that set (we drop the dependence on the process in the notation since it will be clear from 
context), l(-) will be used to denote the indicator function. For any process and a subset A of the 
corresponding state space t(A) will denote the hitting time of that set. When w is a state in the 
space and A is a set of states, we will write {r(A) < K \ uj} to denote the event that, starting from 
a;, the process hits A in at most K steps. 

We will often use the same letter (generally C. D , c or d) for a constant whose value may change 
from line to line. We use x to denote that two quantities that are equal up to universal constants. 
Also O(-), fl(-), ©(•) are used to denote their usual meaning. 


3. Connectivity properties 

Implicit in the statement of Theorem 1.1 is the claim that competitive erosion has a unique sta¬ 
tionary distribution. In this section we formally define the competitive erosion chain and prove this 
claim. 


3.1. Formal definition of competitive erosion. Recall HJ n from Setup 1. We will use HJ n to 

denote both this graph and its set of vertices. For t = 0,1, 2,... let (xi^)<;>o and (yj t+2) ) s > 0 be 
independent simple random walks in HJ n with 


p(4 t] = *i) 


i 




= Z 2 ) 


for all z\,z 2 £ respectively. That is, each walk X® starts uniformly on Ui in and each 

walk Y^ t+ 2 ^ starts uniformly on U 2 n - Given the state S(t) of the competitive erosion chain at time 
t, we build the next state S(t + 1) in two steps as follows. 

S(t+±) = S(t)\J{X® ) } 


where r(t) = inf{s > 0 : 0 5(f)}. Let 


s( t + i) = s( t + l)-{y r ( ; + + |} 


(£_|_ 1 ) 

where r(t + |) = inf{s > 0 : 2 € S(t + |)}. As already stated in (1.9) this formally defines 

the evolution of the process at on Q n U f l' n (see (1.8)). 
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3.2. Blocking sets and transient states. We show that the erosion chain has an unique irre¬ 
ducible class and hence has a well defined stationary measure. 

Definition 2. Call a subset A C U n blocking if HJ n \A is disconnected and the subsets Ui >n \A 
and U2,n\^4 lie in different components. 

Definition 3. For two disjoint blocking subsets A, B C U n we say that A is over B if 

(1) A and Ui, n \-£> lie in different components of U n \B; and 

(2) B and U 2 ,n\A he in different components of U n \A. 

Lemma 3.1. The competitive erosion chain has exactly one irreducible class. Moreover any a E 12 
that has a blue blocking set over a red blocking set is transient. 

Proof. We first claim the following: there exists <r* E 12 n such that for * = 1,2, crf 1 ^) is a connected 
set containing Uj )n . This is true in the continuous setting clearly by taking the sets to be U a and 
(see (1.7)). To find such subsets in U n we look at U a n U n and its complement. These sets might 
not be connected and might not have the exact number of points. However using the smoothness 
of U, and considering the sets bounded by a lattice path approximating the boundary of U a and 
making minor perturbations adding and removing only 0(n ) many vertices gives us such sets. The 
details are omitted. 

Now to prove the first part notice that from any a one can reach cr*. Since in the target configuration 
there is exactly one blue component B and red component R, we look at the closest vertex with a 
value 2 from Ui „ in B. Since there is a blue path from Ui n to that point the Markov chain allows 
us to change it to 1 and similarly at the other end in R. Thus we are done by repeating this. 

To prove the second part we prove that starting from cr* one cannot reach any configuration like 
the one on the right which has one blue blocking set over another red blocking set. We formally 
prove this by contradiction. Let cr be a configuration with a blue blocking set B over a red blocking 
set A. Now the erosion chain evolves as <7o, oq/ 2 > ° 3/2 • • • where for every non-negative integer k 

cr*; E 12 and crk+i /2 £ 12b 

Assume now that there is a path 0o, oq/ 2 , 04 , 03/2 • ■ ■, = a. Let t' and t" be the last times along 

the path such that at least one vertex of A is blue and similarly B contains at least one red vertex 
respectively. If such a time does not exist let us call it —oo. Since do does not have a blue blocking 
set over a red blocking set, max(r', r") > —oo. 

t' must be a half integer (since at half integers there is one more blue particle) and similarly t" 
must be an integer. Thus t' t" . Next we see that we cannot have t' < t". This is because then 
for all times greater than t' , A is entirely red. Thus no blue walker crosses A at any time greater 
than t' . So B must already be blue at t' and stays blue through till t implying that r" < t' . By 
similar argument we cannot have r' > t" . Hence we arrive at a contradiction. □ 


4. Proof of Theorem 1.2 


The goal of this section is to prove Theorem 1.2 modulo some results whose proofs will be deferred 
to later sections. We first restate formally the crucial statistics mentioned in Section 2 to be used 
throughout the article. Recall from (2.2), 


• The discrete Green function : for any x E U n , 



where for * = 1,2, P x (X(t) E Uj iTl ) is the probability that random walk on U n started from 
x , is in U i^ n at time t. c = c(<5) is a centering constant we introduce to ensure that if 6 is 
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small and n is large then the Green function is close to the function 


64 


7T 


log 


(4.2) 


ip(x) — i 
ip(x) + i 

This follows from Theorem 5.1 and Lemma 5.3. c(<5) can be written in terms of the RBM 
heat kernel on HJ. The explicit form of c does not affect the arguments and appears later in 
( 8 . 2 ). 


The weight function: for o £ U fL, 

w(a) = ^2 G n(:J 
xeB 1 

where Bi(a) is dehned in (2.3). 


(4.3) 


4.1. Technical results. We start by defining the following ‘good set’. Consider the following way 
to divide our domain U according to the geodesics defined in (5.8). Let 


... a_i, ao, ai, 02 ,... 


be such that d( 7 ai , la i+1 ) = , for all i = ... — 1 , 0, 1 , 2 , ... and ciq = f3 (recall the definition Ug 

from (1.5)). Thus number of afs is linear in n. Also away from the points X\,X 2 on the boundary 


the function log 


*l>(z)+i 


hence as long as d('y ai , 


(see (5.8)) is a bounded function with bounded non zero derivative and 
xi,X 2 }) > a for some a > 0 (independent of n), 


|oj - Oj_i| = ©(-) (4.4) 

n 

where the constants in the ©(•) notation depends on a and U. Let us define 

U«:=U ai \U ai+1 . (4.5) 


Given e > 0, define, 


12( e ) := Q n r e \ to be the set of all configurations a £ Q n (4.6) 

such that there exists at most en many negative i's such that ILW (IB 2 7 ^ 0, (see (2.3) for definition 
of B 2 ) and at most en many positive i's such that n B\ 7 ^ 0. Thus in words consists of 
conhgurations where at most 2 en many shells U® contain vertices of the ‘wrong’ color. 


Recall the erosion chain at- As mentioned in Section 2 the next result is one of the key technical 
ingredients of the paper. It shows that the weight function has a positive drift if the configuration 
is outside the set 

Theorem 4.1. Given e > 0, there exists a = a(e) > 0 such that for all small enough 5, and all 
n = 2 m large enough, (depending on 5) if a £ then 

E(u;(cji) — w(a) | uo = cr) > a. 

The proof is quite involved and is deferred to Section 6 . Next we state a few properties of the Green 
function (see (4.1)) and the weight function w(-) whose proofs are deferred to subsection 5.3. The 
following result proves a uniform upper bound of the weight function independent of <5. 

Lemma 4.1. There exists a constant D such that for all 5 small enough and n = 2 m > N(5), for 
all a £ O, |rc(fj)| < Dn 2 . 
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Set 

Description 

Defined in 


all the sites out side an e band are the right color 

(1.2) 

A: 

at most 2 en 2 sites have the wrong color 

(1.11) 

r e 

weight function within en 2 of its maximum 

(4.10) 

%) 

at most en different ‘levels’ have sites of the wrong color 

(4.6) 


Table 1. Different kinds of good sets. 


The above statement is a consequence of the following fact: recall from (4.2) that as 6 goes to 0 the 


z)+i 


Green function approaches the function 44 } 0 g 
function divided by n 2 should approximate the integral of t 
The lemma now follows from the fact that the function log 
Let us denote by 

u>i nax := sup w{a). 
cr£f2 


. Now if n is chosen large enough the weight 
le above function over the blue region, 
is integrable. 


j>{z)-i 


ip(z)+i 


(4.7) 


The next lemma shows that the weight function is uniformly close to its maximum value over the 
set A (Table 1) 


Lemma 4.2. For e > 0. There exists constants <5o(e),C( e ) such that for 5 < 5 q and large enough 
n = 2 m > N(5), 

inf w a > Wmax - C n 2 (4.8) 

creAe 

where £ goes to 0 with e. 


As mentioned above the Green function approaches the function 44 log 
by 5 going to 0. Thus in particular it blows up in the limit 
“logarithmic singularity 
to 0. 




^ z ) + i , as n -A oo followed 
The next lemma established the 
of the limit by determining the rate of blow up of the function as 5 goes 


Lemma 4.3. For all small enough S, 

lirnsup sup \G n (z)\ x | log(<5)|. (4.9) 

n=2 m zev n 

We define another ‘good’ set. Given a number e± let, 

r ei = {a £fi: w(a) > w nmx - 2ein 2 }. (4.10) 

Remark 4.1. Given e > 0 there exists 62 < e < e\ such that for all 5 < 5o(e) and n = 2 m > N(S), 

r e2 c Acr £1 . 

Moreover both ei, €2 go to 0 as e goes to 0. Also there exists £3 such that G( £3 ) C A e . 

Thus the above remark strengthens Lemma 4.2. It says that not only does the weight function stay 
close to its maximum on A e , but A e is roughly a level set of the weight function i.e. the actual 
level sets r e ’s can be approximated by the sets A e ’s and vice versa. 

For proofs of Lemmas 4.1, 4.2, 4.3 and Remark 4.1 see subsection 5.3. 
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4.2. Hitting time estimates. The next result proves a uniform bound on the hitting time of the 
set A e . 

Lemma 4.4. Given ei, e 2 > 0 for all 5 < <5o(ei, e 2 ), f/iere exists C, d such that for all n = 2 m > N(5) 
for all (T £ O, 

P (J (r(A 1 n fi (ea) ) > Cn 2 ) < e _dri2 (4.11) 

E ct (t(Ai n fi( £2 ))) < Cn 2 
where 5 is the parameter appearing in Setup 1. 

Now we show that once the process has hit A e it tends to stay “close” to the set for a long time. 
By Remark 4.1 the next result suffices. 

Lemma 4.5. Given e\ > 0 there exists So(ei) such that for all 5 < 5 0 there exists a positive constant 
d = d(e 1 , 5) such that for all large enough n = 2 m > N(S) and any <7 £ 12 

inf P CT (r(12\T 2ei ) > e dn2 ) > 1 - e”** 2 . 
o-er £1 

The above results follow from a general lemma about hitting times for submartingales. The state¬ 
ment has a few parameters and could be difficult to read. However for subsequent applications it 
would allow us to plug in various values for the parameters. 


Let l aft) be a stochastic process taking values in an abstract set V. Also let g : V — > M be a real 
valued function. Let Tt be the filtration generated by the process up to time t and X t := g(ui(t)). 

Lemma 4.6. [15, Lemma 9] Let A±, A 2 ,ai > 0. Suppose 

|< 7 (u;) | < A\ for all ui £ V (4-12) 

|Xt — Xt- 1 | < A 2 for all t. (4-13) 

Also suppose that B C V is such that for any time t 

EpQ-A- 1 IA- 1 ) > ail (w(i - 1) i B) (4.14) 

then 

i. E„,(r(6)) < Moreover P ul {r{B) > T) < exp ^ a2 4 j j or a ll u £ V such that 
g{w ) > A\ — a 2 for some a 2 > 0 and any T such that a 2 — a{T < 0. 


ii. Now consider the special case when B is a level set i.e suppose for some 04 > 2A 2 , B = 
{a; : g(w) > A\ — 04 } and B' = {u : g{w) < A\ — 204} . Then for all u £ B and all T > 


1 


exp i"^ W) 


+ exp - 


32 A^T ) 


^u{t(B') >T')> 1- 
where T' = exp ^ j . 

Using the above result the proofs of Lemmas 4.4 and 4.5 follows easily. 


a\T 2 \ 


Proof of Lemma 4.4. The proof follows from Lemma 4.6 i. The stochastic process we consider 
is the erosion chain a±. Xt = w(&t) is the weight function defined in (4.3). We make the following 
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choice of parameters. 


A 1 = Dn 2 
A 2 = C log(h) 
a\ = a(e) 

02 = 2 Dn 2 , 


appearing in Lemma 4.1 
for a large enough universal constant C 
appearing in Theorem 4.1 
(3D + l)n 2 


a i 

(4.12) is satisfied satisfied by Lemma 4.1. That (4.13) is satisfied by our choice of A 2 follows from 
Lemma 4.3. Thus by Lemma 4.6 i. the result follows. □ 


Proof of Lemma 4.5. The proof follows from Lemma 4.6 ii. Let 

X(t ) = w(a t ) 

where w(-) is the weight function defined in (4.3). We make the following choice of parameters: 

B = T ei ,B' = r 2 ei ,a 4 = ei n 2 ,T = n 2 . 

A-i = w ma . x (see (4.7)), 

A 2 = Clog((5) for a large enough universal constant C , 
ai = 0 ( 63 ) appearing in Theorem 4.1, 

The drift condition (4.14) is satisfied by a\ = a(e 3 ) where a(-) appears in Theorem 4.1 and £3 in 

. / 2 2 \ 2 
min(£pa^)n 

Remark 4.1. Now by the above choice of parameters T’ = re . Thus by Lemma 4.6 ii. for 

all a £ T ei , 


r(r(r 2ei )>r , )>l-[e-T5-+e- 


.Air 


Hence the proof is complete. 


□ 


4.3. Proof of Theorem 1.2. We will use the following general lemma relating hitting times and 
stationary measure. The result roughly says the following: consider any Markov chain and subsets 
A C B of the state space . If the hitting time of the set A is uniformly “small” compared to the 
exit time of B after hitting A, then B has large stationary measure. 


Lemma 4.7. [17, Prop 1.4] Let ut(-) be an irreducible Markov chain on a finite state space A4. 
Suppose A,B C M.. Let 


Then, 


t\ = maxEj,(r(H)) 

x&B 

(4.15) 

t 2 = nrinE x (t(B)) 

xEA 

(4.16) 

7 t(B) < ———. 
v ' _ h + t 2 



where n is the stationary distribution of the Markov chain w(-) on Ai. 


We now finish off the proof of Theorem 1.2. Notice that by the lower containment in Remark 4.1 
it suffices to prove that for given small enough e for all large enough n = 2 m 

vr(T 2e ) > 1 - e-™ 2 


(4.17) 
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for some c = c(e, <5, U) > 0. The proof now follows immediately from Lemma 4.7 with the following 
choices of the parameters: 

A = T e ,B = n\r 2e , ti = din 2 , t 2 = e d2n \ 

where ci, c 2 , d\, d 2 are chosen such that the hypotheses (4.15) and (4.16) of Lemma 4.7 are satisfied. 
Lemmas 4.4 and 4.5 allow us to do that. □ 

5. Technical Preliminaries 

In this section we will develop some preliminary tools needed for the proof of Theorem 4.1. 

5.1. Discrete Green Function. Recall that we consider continuous time random walk on U n 
(Section 2.3 ii). Call it X(t). For x, y G U n let 

P x(X(t) = y) (5.1) 

denote the chance that the random walk on HJ n starting from x is at y at time t. For notational 
simplicity we suppress the n dependence in P since the graph will be clear from context. Similarly 
for any set A C U n , let P x (X(t) G A) denote the chance that the random walk is in A at time t 
starting from x. Recall the definition of the function G n on U n from (4.1): for any x G U n 

2 n 2 r°° 

Gn(x) = / \P x (X(t) G Ui >B ) - F x (X(t) G U 2 ,n)}dt - c. (5.2) 

|Ul,re| Jo 

where c = c(S) is some constant explicitly mentioned in (8.2). The fact that the integral in G n is 
absolutely integrable follows immediately from the following lemma. 

Lemma 5.1. [ , Lemma 3.1] Given U as in Setup 1 there exists a constant D = D(U) and a time 

T = T(U) such that, for all large enough n, 

sup \F x (X(t) G U lin ) - P*(X(t) G U 2iri )| < 2e~ m 

x£U n 

for all t >T. 

The above is a straightforward consequence of the mixing rate of random walk and the fact that 
since by choice |Ui in | = |U 2 , n |, (see subsection 1.5 ii.), 

Tmy(Ul,n) = TRVu(U2,n) ; 

where itrw is the stationary measure for the random walk on U n . 

For any function / : U n —> M define the laplacian A/ : HJ n — > M: for any x G HJ n 

A/(x) = f{x) - f(y) ( 5 - 3 ) 

x y^x 

where d x is the degree of the vertex x and y ~ x denotes that y is a neighbor of x. 

Lemma 5.2. Consider the function G n (-) on Un. Then, 

A(Gn) = —j-( 1 (U 1 ,n) - 1 (U 2;77 .)), 

where for any subset A C U n , 1(^4) denotes the indicator of the set A. 


(5.4) 
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Proof. Proof follows from definition of G n and looking at the first step of random walk which by 
definition is of expected duration 2 f 2 . Thus we have 

G n (x) = W^(l(U 1>rl ) - l(U 2 ,n))(x) + j- £ °n(y) 

' ,n ' X yr^X 

and hence the lemma. □ 

5.2. Conformal invariance of Green function. In this section we state a key technical result 
establishing conformal invariance of the Green function. We define the following function on the 
smooth domain U : 

/; = aleak) ( 1 d , 2 )-10 , l)). (M) 

where HJi and U 2 are defined in Setup 1. Again 16 is a constant that appears for the same reason 
that 64 appears in (1.5) and is not important. 


Recall the functions f and if from (1.4). Now let us call ififJi) = Aj. Then 

16 


f°<t> = 


area(Ui) 

For any U as in Setup 1, define the function G* : U 
that (f(y) = z, then 


(l(A 2 ) — l(Ai)). (5.6) 

M such that for all z G U if y G D is such 


G,(z) = - 

7T 


ICI<1 


/ o 0(C)I0 , (C)I 2 iog(|(C - y){ 1 - Cy)\ 2 )d£dri, 


(5.7) 


where £ = £ + if]- The above expression is well defined on U since by Section 1.5 i., the maps 
(j),^ can be extended to neighborhoods containing B,U respectively. Notice the dependence of G* 
on 6 through /. However for brevity we choose to suppress the dependence on 6 in the notation. 
The definition of the function might be a little mysterious at this point (the reason for our choice 
is explained in Section 8.1). Nevertheless Lemma 5.3 shows that as <5 goes to 0, the function 
approaches the function log , up to a multiplicative constant. Since the above function 

is conformally invariant by definition, this establishes that the function G*(-) is “asymptotically” 
conformally invariant. 

Even though the function G n in (5.2) is defined on the graph U n , in this section we use interpolation 
to think of it as a function on the closure of the whole domain, HJ. We are now ready to state one of 
the main convergence results of this paper establishing asymptotic conformal invariance of G n (-). 
We start with the interpolation scheme. We follow the scheme mentioned in [11, Pf of Theorem 
2.12] ( also appears in [13, Pf of Theorem 2.2.8]). 

For all x G ^Z 2 \U n , define G n (x) = 0. Now having defined G n (x) for all x, € we define it on 
C and hence by restricting on U. This is done by interpolating G n {x ) by a sequence of harmonic 
extensions along simplices, 

i. First extend it along the edges using the value on the verices in -Z 2 . 

ii. Then extend it to the squares using the value on the edges. 

Thus the function is extended to the entire complex plane C. By abuse of notation in the following 
result we call the extended function G n {-) as well. 


Theorem 5.1. [16, Theorem 3.1] For all small enough 5 , 

lim sup | G n (z) — GJz)\ = 0. 

m—> 00 _ 1 v 7 \ \ 

n=2 rn ZE U 
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The basic outline is sketched in Section 8 . For the complete proof see [ ]. We end this section 

by making a few observations about the function G*(-) and hence by the aforementioned theorem, 
about the asysmptotics of the function G n . 

Let us define the following one parameter family of curves for any U which form the boundary of 
the sets defined in (1.5). For any f3 G M define 

if(z) — i 


64 

7/9(U) = {z E U : — log 
7r 


ip{z) + i 


= (3}. 


(5.8) 


We will often suppress the dependence on U and denote the above just by 7 ^ for brevity, if the 
underlying domain is clear from context. It is well known that for any /3 E M, 7 ^ is a hyperbolic 
geodesic. In case of the disc, 7 ^’s are circular arcs orthogonal to c© symmetric with respect to — i 
and i. For more details regarding the hyperbolic metric see [1]. 

The following is a uniform convergence result of G*(-) away from the points xi,X 2 - 
Lemma 5.3. Let 4J be as in Setup 1. Given a < 1, 

ip(z) — i 


lim 

<5—5-0 


sup 

z6®: 

d(z,xi)>S a 

d(z,X2)>S a 


G*(z) - — log 
7 r 


ip(z) + i 


= 0 


Proof. Recall that the map <j) is bi-Lipschitz by subsection 1.5 i. Thus it suffices to show that 

= 0. 


lim 

5-s-O 


sup 

G* 0 <f)(z) - — log 

z — i 

z + i 

,-i)>5 a 

7 r 


Recall the sets Ai,A 2 C 


and hence 


d(z,i)>S a 

from (5.6). Note that by hypothesis 
Ui C B(xi,2<5)nU 


Ali C B(-i,C6)nB (5.9) 

for some constant C ( where C /2 is the Lipschitz constant of ip). Similar containment holds for U 2 
and A ‘2 ■ Now notice that by the change of variable formula for i = 1,2, 

1 


area(Ui) J A . 

These facts along with (5.7) clearly imply that 

64 

G* o (f(z) = — log 

7T 


fomWiQfdfdp =16. 

\z — i\ + 0(5) 

\z + i\ + 0(5) 


(5.10) 


(5.11) 


where the constant in the O(-) term depends only on O in (5.9). Thus 


64 

G* o cf>(z) -log 

7 r 


z — 1 


z + i 


= O 


log ( 


1 + 


0(5) 


z — i\ 


+ 


log 1 + 


I z + i I 


Since a < 1 we have lim — = 0. Thus sending 5 to 0 we see that the RHS uniformly converges to 
5-s-o o a 


0 on the set (zGD: min(|z — i\,\z + i\) > 5 a }. Hence we are done. 


□ 
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Since the function log 


z+i 


, has logarithmic singularities at i and — i as 5 —> 0 by Lemma 5.3 the 


function G*(z) blows up near the singularities x\,X 2 - The following lemma determines the rate at 
which the blow up occurs. 

Lemma 5.4. For a < 1 and all small enough 5, for all z G HJ such that d(z,x i) < 5 a 

G*(z) ~ |log(5)| 

and similarly if d(z,x 2 ) < 5 a , 

G*(z) x — | log(<5)|. 

The constants in the x notation depend on a only. 


Proof. We only show the first case since the argument for the other case is similar. Hence we 
assume d{x 1 , z) < 5 a . Also using the bi-Lipschitz property of f it suffices to show that for all z G D 
such that d(z, — 1 ) < 5 a , 

{G* o 4>)(z) x | log(<5)|. 

Using (5.7) we do it in two steps. We first show the following: 

W- 

-L-i / ( 1 (A 2 ) - l(.4i))(C)|0'(C)| 2 log|C - z\d( - (I log(J)l). 

area(Ui) ir J\q<i 

Clearly 

-Lpri / 1O2)(C)I0'(C)| 2 l0g|{ - 2|<i< = 0(1) 

area(Ui) 7r 

since d(z,A 2 ) > d(xi,cc 2 )/2 for 5 small enough. Since if is bi-Lipschitz, A\ is contained in 
a ball of radius CS and contains a ball of radius c5 for some constants C, c not depending 
on <5 which implies that area(Ai) = 0(<5 2 ). Also \cj)'\ is bounded away from 0 and 00 . 


Thus for any z£Dwe have 

1 1 f 


area(HJi) 7 r 


i(^4 1 )(C)|^ / (C)l 2 iog |C - 


= o 


= o 


1 

52 

1 


rCS r 

'0 J\g-z\=r 

r C5 


log |C — z\dddr 


I r \og(r)dr 


1 0 


= O (| log(<5) |). 

The lower bound follows immediately since d(z,A\) < 5 a . 


(ii). We now prove a similar bound for the term 

1 1 f 


area(Ui) 7r J\q k1 
Again it is easy to check that 

1 1 r 


area(HJi) 7r J | C |<i 
T o find the exact order of the term 

1 1 f 


area(Ui) 7r 


1(^40 - i(A 2 ))(OI0 , (C)l 2 iog|i - C^MC- 
l(A 2 )(C)|^(C)| 2 log|l-^|dC = 0(l). 
i(A 1 )(C)|0 , (C)l 2 iog|i - 
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recall from Setup 1 d(Uj,<9U) = 0(5) and hence d(Ai,dU>) = 0(5). Thus for any f £ A\, we 
have C5 < |1 — C, z \- Also since d(z,x i) < 5“ by hypothesis jl — C,z\ < 0(5 a ). Therefore, 

—L-l / (i(Ai))(0l/(0l 2 iog|i-£K-(|logWI). 

area(Ui) 7 T y| C | <1 v 

Thus putting (i) and (ii) together we are done. □ 

5.3. Proofs of some earlier statements. In this section we provide the proof of some of the 
statements stated earlier. 


Proof of Lemma f.3. It follows immediately from Lemmas 5.3, 5 


.4 and Theorem 5.1. 


□ 


Proof of Lemma f.l. Clearly by (5.7) and Lemmas 5.3, 5.4, G*(z) is a bounded continuous function 
on U (where the bound is 0(| log(5)|). Thus approximating integral by Riemann sums we get 

\G*(z)\dfdr] = lim W \\G*(z n )\ 

n—>• oo z —* Ti 
^nGUn 



where z = £ + irj. The proof will then immediately follow from Theorem 5.1 if we can show that 


lim sup / |G*(z)| df,dr] < oo. 

<5 —>-0 J u 


Now notice that 


lim [ \G*(z)\dfdi] = f — 


<5—>0 


'U 


JV 71 
64 r 

77 Jo 

< oo. 


log I 


ip(z) — i 


ip(z) + i 
\(tf{z)\ 2 log| 


dfdr] 


(5.12) 


z — i 


z + i 


dfdr] 


The first equality is a consequence of Lemmas 5.3 and 5.4. The second equality is just the change 
of measure formula. The last inequality holds since |<^ / (z)| is bounded and the function log| 'jfff | is 
an integrable function on the disc. Thus we are done. □ 


Proof of Lemma f.2. Let us denote by a* a configuration where every vertex in U( a ) except possibly 
0(n) many, is colored blue. The 0(n) correction is needed since every configuration has exactly 
|_a|U n |J many blue vertices. By Lemma 5.3 for all z £ HJ such that min(d(,z, aq), d(z, X 2 )) > 6 a we 
have, 


64 


G*(z) = — log| 

7 r 1 yj{z 


if{z) — i | 


+ i 


+ o(l)- 


(5.13) 


By Lemma 4.3 on the remaining set of measure 0(5 a ), one has |G*| = 0{\ log(5)|). This along with 
Theorem 5.1 implies that given any number c for small enough 5 and n = 2 m > N{5) we have, 


»max - n 2 — f log| | dfdri 

7T Jv a 


1 if(z) + i 
|iu(<7*) — W 


max ^ 


cn 


cn 


(5.14) 


From Lemmas 5.3 and 5.4 it follows that the function G*(z) (as 5 goes to 0) is uniformly integrable 
i.e. given any a there exists b = 6 ( 0 ) (b(a) —> 0 as a goes to 0 ) such that, 
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where the supremum is taken over all subsets A of U of measure at most a. Using this and Theorem 
5.1 we thus have the following : Given any constant a there exists a constant b = b(a) such that 
for all small enough S and n = 2 m > N(5): 





where b = b(a) tends to 0 as a tends to 0. Thus for any a e A e by (5.15), 

| w(a) — w(a *)| < b(2e)n 2 . 

Hence we are done by (5.14). 


□ 


Proof of Remark f.l. The upper containment follows from (4.8). The lower containment fol¬ 


lows from (5.13) and Theorem 5.1. 


□ 


We are now ready to prove the main technical result of this paper. 


6. Proof of Theorem 4.1 


The proof is quite long and involved and constitutes the technical core of the paper. The argument 
has to be divided into cases depending on a. Thus we split the result into two theorems. We will 
suppress the a dependence in the notations we define subsequently for brevity whenever there is 
no scope for confusion. We start with some definitions and notation. 

Definition 4. Given tr 6 H U fl' (see (1.8) for definition) for i = 1,2, let Ri = Ri(cr) be the set of 
all points in 4J n of color i reachable by a monochromatic path of color i from a point in U j n ( we 
do not allow the endpoint to be the opposite color). Note that Ri and i ?2 are disjoint. 



Figure 6. The blue connected region containing the blue blob is R± and similarly 
i? 2 - Note that the two red islands and the blue island are not included in either set. 
Also recall (even though not in figure) that if a part of the blue blob is colored red, 
it is not included in the set R,\. 

Clearly by definition 


E H<u) - w{ctq) \<t = (t 0 ]= E[G n (X r{Rl{ao)) ) - GnQ4(R 2( <x 1/2)) )] (6-1) 


where X r ( Rl ( ao \\ is the point at which the random walk exits Ri(cro) and similarly the second term. 
The expectation in the first term is over /ii (starting distribution of the blue random walk, see 
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Setup 1.) Note the expectation in the second term is over fi 2 as well as the random intermediate 
configuration oy/ 2 , (see (1.10)). 

For any x £ U, a > 0 define 

C(x , a) := dB(x, a) n U, i.e. the part of the boundary inside U. (6-2) 

Recall the weight function w(-) from (4.3). We now state the first theorem towards the proof of 
Theorem 4.1. The theorem roughly says that if either diameter of R\ is “small” or R 2 connects a 
region “close” to X 2 to a region “close” to x\ then the drift in Theorem 4.1 is like | log(<5)|. See Fig 
7 and 8. 

Theorem 6.1. Let U be as in Setup 1. Then there exists constants 

0 < 02 < Oi < 1 

and 5o such that for 5 < do and all n = 2 m > N{5) if a G Cl is such that either 

i. R,\ does not intersect C(yi,S 2 ) or, 

ii. R\ intersects C(yi, 62 ) and R 2 intersects C(yi, 5\) or 

switching the roles of R± and R 2 in the above two cases 

Hi. R 2 does not intersect C(y 2 ,d 2 ) or, 

iv. R 2 intersects C(^ 2 , 82 ) and R\ intersects C(y 2 ,S\) 


then 

E[u;(cji) - w(ao) \ cr 0 = a] x | log(<5)|. (6.3) 

The reason for the above is roughly the following: (we only discuss the first two cases.) 

• In case i. diameter of R\ is small implies that the blue random walk stops before exiting a 
small ball. Thus the first term in (6.1) is roughly | log(<5)| by Lemma 5.4. Now by standard 
random walk estimates one can show that the red walk is likely to exit R 2 before coming 
close to x\. Thus the second term in (6.1) is much smaller and hence the result follows. 

• In case ii. by hypothesis there is a red path which connects \$ 2 ,n to a small ball around y\ of 
radius 8 a for some a. Then it becomes likely that the blue random walk starting from Ui iTl 
will hit this path and hence exit R± before going too far from x\ again makes the first term 
in (6.1) roughly |log(5)|. The rest of the arguments are then the same as in the previous 
case which shows that the second term in (6.1) is much smaller. 

Before providing the formal proof we first state a lemma regarding random walk on U n . Let P x (-) 
denote the random walk probability measure started from x. 


Lemma 6.1. [16, Lemma 5.2] For all U, x\ as in Setup 1, 


lim lim sup sup 


e— *0 rwoo creQUO' LzeU„\B(a:i,|) 
where B 1 = B\(a) was defined in (2.3). 


sup ¥ z {T(B(xi,e)) < t(B 1 )} 


= 0 


The above lemma says that uniformly from any point z at distance 1/2 (any constant would work) 
from x\ and all configurations a G Cl U fL the chance that the random walk does not hit B\ before 
reaching at distance e from x\ goes to 0 as e goes to 0. Note that since B\ is the set of all blue 
particles by definition ttrw(B 1 ) = a — o(l) where the o(l) goes to 0 as n goes to infinity. The proof 
of the above lemma for general sets of measure at least a appears in [ ]. We are now ready to 
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X\ 


Figure 7. Illustrating the case when the diameter of R\ is less than 6 2 and hence 
the blue walker stops within B(x 1 , im¬ 


prove Theorem 6.1. 

Proof of Theorem 6.1. 

By symmetry it suffices to prove the theorem for the first two cases. 

Proof of case i. Let 03 < 02 < a\ < 1 be constants to be specified later. Let 63 = <5 “ 3 and 
p(<5 3 ) be the probability that the random walk started uniformly from U 2 jn hits I?(xi,< 5 3 ) before 
hitting B\(a) (the set of blue sites) which has size aarea(U)n 2 + O(n). Now observe that in this 
case there exists universal D\, D 2 > 0 (not depending on S, n) such that for all small enough 5, and 
large enough n 

E(u;(cti) - w(cjo)|cjo) > Z>i|log(<5 2 )| - [p(5 3 )D 2 \ log(<5)| + (1 - p(5 3 ))D 2 \ log(<5 3 )|]. (6.4) 

To see (6.4) we first recall (6.1). 

E['in(or) - w(a 0 ) \ a = <r 0 ] = E[G n (X r{Rl ^ o)) ) - C? n (X r(Ba((Tl/2)) )]. 

By hypothesis R\ does not intersect B c (xi,5 2 ). See Fig 7. Hence the blue random walk stops 
before exiting B(xi,5 2 ). Thus from Theorem 5.1, Lemma 5.4 it follows that the first term is at 
least D i| log(< 52 )|- Now the second term is at most 0(\ log(<5)|) by Theorem 5.1 and Lemma 4.3 if the 
red random walk enters the ball i?(a?i, 5 3 ) before hitting - 61 ( 07 / 2 ) which happens with probability 
at most p(d 3 ). Otherwise it is at most D 2 \ log(<5 3 )| by Lemma 5.3 and Theorem 5.1. 

Since p(<5 3 ) goes to 0 as 6 goes to 0 we can suitably choose a 3 < a 2 to be done. Note that we do 
not need 01 for this argument. It appears in the proof of the next case. 

Proof of case ii. By hypothesis there exists a path 7 of red vertices connecting B(xi,5\) and 
B c (xi, S 2 ). Thus it should be likely that the blue random walk started from Ui n will hit 7 before 
exiting B{x 1 , S 2 ). This is precisely the content of the next lemma which roughly says if a connected 
set A of large enough diameter is close enough to Ui n then random walk starting from Ui n is more 
likely to hit the set A before exiting a large enough ball. 

Lemma 6.2. [16, Lemma 5.5] Let 0 < e\ < e 2 . Assume A C U n is a connected set such that 

d(Ui in ,H) < ei. 
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X\ 


Figure 8. Illustrating the case when f? 2 connects B(y 2 , 5) to a point in the neigh¬ 
borhood of X\. 

Also assume An (\J n \B(yi, 62 )) k 0. Then 

sup F > x(r(U n \S(yi,e 2 )) < r(A)) < C Iog ( e i), 
x£Ul, n 

for some C = C(U) < 1 independent of n. 

Thus by the above lemma, 

inf P x (r( 7 ) < t(B c (xi, 52 ))) > 1 — C' log ^ 2// ' 5l ' ) (6.5) 

xeUi,„ 

for some C = C(U) < 1. Thus 

E(w(oi)-w(o- 0 )ko) > Di\ log(<5 2 )| -C' log( ' 52/<5l) T> 2 | log(5)| - [p(<5 3 )| log(<5)| + (1 -p( 6 3 ))D 2 \ log(d 3 )|] • 

( 6 . 6 ) 

This is because G n (X T ( Rl r a0 ))) is at least |log(d 2 )| if the blue random walk stops before exiting 
S(xi, 5 2 )) which happens with chance at least 1 — C 1 °g(‘W' 51 ) by (6.5). Otherwise it is at least 
—0(| log(<5)|) since that is the minimum value G n takes by Theorem 5.1 and Lemma 4.3. Bound 
on G n (X r ( K2C - i/2 )) is exactly as in case i. 

Now one can suitably choose constants a 3 < a 2 < a± such that the RHS in both (6.4) and (6.6) are 

0(| log(<5)|). Thus like case i. the proof follows by suitably choosing aj’s for i = 1,2, 3. The proof of 

Theorem 6.1 is hence complete. □ 

Now to finish the proof of Theorem 4.1 we have to consider the cases not considered in Theo¬ 
rem 6.1. Thus we first make the assumption, 

Assumption 1. : a does not fall in any of the four cases in the statement of Theorem 6.1. 

Throughout the sequel till the completion of the proof of Theorem 4.1 we will work 
under the above assumption. 

That is both f?i,f? 2 extend beyond B(x\, 62 ), B(x 2 , <5 2 ) respectively and also f? 2 does not reach 
B(xi,di) and similarly R± does not reach B(x 2 ,Si). Thus under Assumption 1, 

B(x ll 6 1 )f]R 2 = 0 (6.7) 

B(x 2 , 6 1 )f]R 1 = 


0 . 
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Notice that for i = 1,2, 


C B(yi, 5\) 

(see Setup 1 for definition of ILJ n ) for small enough 5, since <5i 


8 a for some a < 1. 


Remark 6.1. We now claim the following. Given any small enough 5 for all n > N{$) there exists 
a set A\ iTI C U n such that Ai jTI is connected and, 

B(y 1 , |)nu„c A 1>n c B{ yi , 5i) n U„. 

Similarly one has A 2 tU corresponding to 1/2 ■ This is a simple consequence of the fact that the locally 
near the boundary HJ looks like a half plane, see (9.1). We omit the details. 


Now since 

r 

U i t n C B(l)i, —) Pi U n C A{^ n , 

it immediately follows from definition of Ri that for i = 1,2 

R-i A l n 

are two disjoint connected subsets of HJ n . For technical reasons sometimes it will be convenient to 
work with the following larger set instead of R\. Let HJ n \{i?i J A\ tn } = Ui=i Ci,» where C\.i are the 
connected components of U n \{l?i |J A\ )n )}. Let {R 2 IJ ^.n} C C\,k- Define 

Ri = Ri U Ai iTl Cij. ( 6 . 8 ) 

i^k 

In words R\ is the set R\ J Ai >n along with all the red holes filled in. Similarly define 1?2- Note that 
for i = 1,2, Ri is connected. Note that Ri depend on o which we choose to suppress for notational 
cleanliness. See Fig. 9. 



Figure 9. Illustrating the difference between the sets Ri and Ri for i = 1,2. Note 
that in the latter the holes are filled in. 

Given d C U n , we define two types of boundary of A, 

d ou tA := {y E A c : 3 x E A such that x ~ y} (6-9) 

di n A := {x € A : 3 y G A c such that x ~ y}. 

Definition 5. Let HJ* denote the graph U n along with all the diagonals of the squares that are 
entirely in U n . We will call connected subsets of U* as *— connected subsets of U n . 
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Remark 6.2. Since both R\ and Rf are connected by Remark 6.1 it follows by [30, Lemma 2] that 
both d ou tR\ and di n R\ are *—connected. 

We now proceed towards completing the proof of Theorem 4.1. The proof will be based on an 
electrical resistance argument. We first review some basic facts about energy of flows on graphs. 

6.1. Energy and flows on finite graphs. For a graph G = (V,E), let w ~ v signify that w 
is a neighbor of v. Also let E be the set of directed edges where each edge in E corresponds to 
two directed edges in E, one in each direction (except for self-loops, which correspond to just one 
directed edge in E). A flow is an antisymmetric function / : E —> M (i.e. a function satisfying 
f(w,v) = —f(v,w)). Note that by definition the value of a flow on a self loop is 0. Define the 
energy of a flow / by, 

£ (f) = l f( v ’ w ) 2 - ( 6 - 10 ) 

(v,w)ee 

For any flow / : E —> M define the divergence, div / : V —> M by, 

div/(u) = Y f(w,v). (6.11) 

W~V 

Note that for any flow, 

^div/= Y f( x ’V) + = 0, (6-12) 

x£V x,yev 

■yrux 

since / is antisymmetric. For disjoint subsets A, B C V and a flow / we say that the flow is from A 
to B if di v f(z) = 0 for all vertices except vertices of A and B and the divergences across vertices 
of A and B are non negative and non positive respectively. For more about flows see [26, Chapter 
9]. For any function F : V —> M define the gradient X7F : E —> M by 

VF(v,w ) = F(w) — F(v). 

Recall the definition of laplacian from (5.3). Thus clearly for any F : V —> M, 

[div(VF)](u) = d v .(AF)(v), (6.13) 

where d v is the degree of the vertex v. The next result is a standard summation by parts formula. 

Lemma 6.3. For any function F : V —> M, 

£(VF) = Y F(v)d v .AF{v). 
vev 

The proof follows by definition and expanding the terms. We omit the details. 

We now discuss some energy interpretations of the Green function G n (-) defined on U n in (5.2). 
The next lemma is a standard result adapted to our setting. 

Lemma 6.4 (Thomson’s principle). £(VG n ) < £(9) for all flows 0 on U n such that 

div 9= i=^ 1 [l(U lin )-l(U 2 ,„)]. 

Proof. First observe that by (5.4) and (6.13), the flow VG n has the same divergence condition 
as above. Note the 4 appears since every vertex in Ui, n and IL^n has degree 4. The proof now 
follows by standard arguments. See [26, Theorem 9.10]. We sketch the main steps. One begins by 
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observing that the flow XG n satisfies the cycle law i.e. sum of the flow along any cycle is 0. To see 
this notice that for any cycle 


Xl,X 2 , ...Xk=X 1 


where xfs E U n , 

k— 1 k— 1 

^ G n (xjiXj-\- 1 ) — G n (xi)) — 0 . 

i =1 i=l 

The proof is then completed by first showing that the flow with the minimum energy must satisfy 
the cycle law, followed by showing that there is a unique flow satisfying the given divergence 
conditions and the cycle law which implies that VG n is the unique minimize! - . 

□ 


= A 


For notational cleanliness we introduce the following notation. 

1 

|Um 

1 

/ 

Thus rewriting (5.4) in this notation we get 

A G n = f. 


|l(Ui, n ) 

rl(U 2 ,n) 


h 

A - h- 


(6-14) 

(6.15) 

(6.16) 

(6.17) 


6.1.1. Stopped Green functions. For i = 1,2, and a E D U fl' define 

G i: n(x ) := Gi, n (cr)(x) := 2n 2 E x [ f(X(t))dt, (6.18) 

J o 

where X(t ) is the continuous time random walk on 4J n , as defined in (5.1) and 

n ■■= n{a) := (6.19) 

i.e. hitting time of i7?( a) for the random walk Xt, where R = Ri(a) was defined in Definition 4. 
Similar to (4.1) we again suppress the dependence on 5 in the notation Gi tU . Now notice, 

A G i} n(<j)(x) = f(x) for x E Ri(a) (6.20) 

G ! j,n(o')lu ri \Ri(cr) = 0. (6.21) 

Recalling notation for the energy of a flow from (6.10), we state the following key lemma, which 
translates Theorem 4.1 to a statement about difference in energies of flows. 


Lemma 6.5. Under Assumption 1, 

E(w(cji) - w(a 0 ) | cro = cr) > |[ £(VG n (a 0 ) - £(VG i,„(o - 0 )) - £(VG 2 , n (0o))]- (6.22) 

Thus using the above result our approach to proving Theorem 4.1 involves lower bounding the RHS 
in the above equation. 

Before proving Lemma 6.5 we make a few more observations about the functions Gi }U to be used 
later. First notice that if X^ was a discrete time random walk then (6.18) is exactly the same as 

T» —1 

E* J] f(X k ), 

k =0 


(6.23) 
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where r,; is now the exit time of the discrete random walk. 

Thus for the rest of the section for notational convenience we will switch to the sum notation and 
hence think of the random walk as discrete time. 

The next remark gives us a crude upper bound for the stopped Green functions Gj >n . 

Remark 6.3. For i = 1,2 and any dGilUll', 

limsup sup Gi >n (x ) = O 

n—2 rn xeVi :n 
m—too 

A sharper upper bound of 0(log(|)) can be obtained with a little more work. However for our 
purposes this will suffice. 



Proof. Recall Rj(cr) from (2.3). Since Ri C Bi it follows that t(R£) < t(B£). As a simple conse¬ 
quence of the fact that the Loo mixing time of the discrete time random walk on Un is 0(n 2 ) (for 
a proof see [16, Lemma 5.4]) and the fact that both B\,B 2 have @(n 2 ) vertices, one gets that for 


* = 1 , 2 , 


sup E(r(Rj)) = 0(n 2 ), 
zeu„ 


Since B\ = B£ from the above discussion we get that, 


sup E(t;) = 0(n 2 ), 

a:gU n 


where r* appears in (6.19). Now the remark follows once we observe that by (6.23), 

G i,n (x) < rE x (Ti). 


□ 


Remark 6.4. By (6.7) under Assumption 1, R\ does not intersect U 2 ,n which is the support of f 2 
and similarly R 2 does not intersect Ui iTl which is the support of f\. 


Thus (6.23) becomes 


Gi, n (a)(x) 

G2, n (<r)(x) 


and (6.20) and (6.21) become 


n ( o -)-1 

E,[ fi(X k )] 

k =0 

T 2 (< t )-1 
k =0 


AG i>n (a)(x) = 

G i,n{ a )\v n \Ri(a) = 
An easy but useful observation is the following: 


fi(x ) for x e Ri{a) 

0. 


(6.24) 


(6.25) 

(6.26) 


Lemma 6.6. Under Assumption 1 for any <r£f l, 

L[(G 2 ,n(cri/ 2 )(x) I £J 0 = a)] < G 2 , n (cr)(x ) (6.27) 

for all x G HJ n where the expectation in the first term is over the random intermediate state oq/ 2 , 
(see (1.10);. 
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Proof. The proof immediately follows from (6.24) since 72 ( 01 / 2 ) < 72 ( 00 ) as 0 1 / 2 has one more blue 
particle than <tq- □ 


6.1.2. Proof of Lemma 6.5. Let us denote the random walk starting uniformly from Uj n by Xk 
and the random walk starting uniformly from U2 n by Y&. Recall (6.1), 

E(w(ai) - w{a 0 ) \ 00 = 0 ) = E(G n (X Tl{ao) ) - G n {Y T2 ^ 1/2) ), 

where the expectation on the right is over the random walks A&, Y^. We start with the following 
telescopic sum, 


GniX-nfo)) ~ Gn(Xo)+[G n (Xi)-G n (Xo)]-\-[G n (X2)—G n (Xi)]+.. .+[G r) .(X ri ( CTo ))—G n (Y ri ( (ro )_ 1 )]. 
Now clearly 

E(G n (A i+1 ) - G n (Xi) | Xi) = -AG n (X i )l(r 1 (0 O ) > *) 

where X% is the filtration generated by 00 and the random walk X up to time i. Taking expectation 
on both sides we get, 

E(G„(A Tl(fTo) )) = E[G n (A 0 ) - Y AGn^w).!)]. 

fc> 0 

A similar equation holds for E (G n (X V2 ^ rTi/2 ^)). Plugging in both of them we get, 

E(G n (A Tl((7o) ) - G n (Y T2{ai/2) )) = E[G n (A 0 ) - Y AG' n (A {feAri((To)) _ 1 ) - G n (Y 0 ) + Y AG n (Y (fcAT2((r 

fc>0 fc>0 

(6.28) 

Now 

E(G n (A 0 ) - G n (Y 0 )) = Y G n (x)(fi(x) - f 2 (x)) = (G n , A G n ) = \s(VG n ). (6.29) 

xeVn 

where all inequalities but the last are by definition. The last inequality is by Lemma 6.3 and the 
fact that all vertices in Uq n U U 2 , n have degree 4. Also, 

E.x ^ ^ AG n (X( fcAri ('o- n ))_l ) E x ^ ^ fl {-X(kATi (q~o))_l) Gi } n(ao)(x), 

k> 0 k> 0 

where the first equality is by (6.25) and the second equality is by definition (6.24). Similarly 

AG ri (Y(fc Ar2 ( 0 - 1 / 2 ))_i) = — E(G2,n(°‘l/2)(y))- 

fc >0 

Hence by the above discussion, notice in (6.28), 


eEac»(x ( b „ wh )] = e eE AG n (A( fcAri ( cro ))_ 1 ) | Ao] 

k >0 k>0 

= E[Gi, n (ao)(A 0 )]. 
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Similarly we have EE fc > 0 AGnC^mj-i)] = -E[G 2 , n (cr 1/2 )(:Fo)]- Thus 


E(to(fJi) - w(a 0 ) | ct 0 ) = 


> 


XI G n(x)fl(x) - XI G n( X )M X ) ~ XI G ^A a o)( x )fl( x ) - G2 1 n{cr 1 /2)(x)f 2 (x)] 

Ul,n U 2 ,n Ul i71 U2,n 


XI G n( x )fl( x ) ~ XI G n( X )h( X ) ~ XI G l,n(°o)( x )fl( x ) ~ XI G lA a o)( x )M x ) 

Ul,n U2,n Ul,n U2,n 


where for the inequality we replace the last term using Lemma 6.6. Now by Lemma 6.3 and (6.25) 

Y, Gl,n(x)/l(x) = ^(VGl,„) 

XI Gf 2 , n (s)/ 2 (s) = ^(VG 2 , n ). 

3 JGU 2 n 


Thus we are done using these and (6.29) in the above inequality. 


□ 


We need a few more definitions and results to use the inequality in Lemma 6.5 to prove Theorem 
4.1. Recall the standard gluing operation on any multigraph G = ( V,E ) where certain subsets of 
vertices are identified (“glued”) and they act as a single vertex. Any edge between two identified 
vertices now act as a self loop in the glued graph. For more on glued graphs see [27]. Also in the 
sequel we will often specify a graph by just referring to the set of vertices which will be a subset of 
the vertex set of U n . The graph then would be the sub graph induced on the set of vertices from 

U„. 

For the following definitions, recall (6.8) and (6.9). 

Definition 6. i. For i = 1,2, let R° ut denote the graph with vertex set Ri U d ou tRi with the 
following subset of vertices, 

|Ri U d out Ri} \Ri 

glued. 

ii. For z = 1,2, let R denote the graph with vertex set Ri with the following subset of vertices, 

{Ri\Ri}[jdi n Ri 

glued. 

The set of glued vertices now acting as a single vertex will be denoted by 0 (the underlying graph 
will be clear from context). 

With the above definition, V(Gj in ) is a flow on R° ut . (6.25) and (6.26) are the same as 

A(Gi tn )(x) = fi for all x £ Ri (6.30) 

Gi, n (t>) = 0. (6.31) 

Analogous to (see (6.24)) for i = 1,2, we define G* n to be the function on Rf 1 such that 

OO 

G ln(x) =® X J2 A X k) l((fc < r(Ri U d in Ri))), 
k =0 


(6.32) 
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i.e. the random walk is stopped on hitting or the boundary di n Ri. It is easy to verify that 

A(G* n )( x ) = fc for all x £ Ri\d in Ri, (6.33) 

G*» = 0. (6.34) 

Notice that by (6.12), the above information determines the laplacian (equivalently divergence) at 
0 in the glued graphs R° ut , R]'f 1 respectively. 

Let us form the following abbreviations 


£(VG n ) = £. (6.35) 

Also for i = 1.2, 

f(VG iin ) = £i (6.36) 

£(VG* J = £*. (6.37) 

For brevity we write (6.22) in terms of the above abbreviations, 

E(tc(fJi) - w(a 0 ) | cr 0 ) > -[£ - £i - £2]- (6.38) 

Remark 6.5. Notice that the arguments in the proof of Lemma 6.4 are not special to the potential 
function G n , i.e. [26, Theorem 9.10] also gives us the following two facts: 

i. £i = min£{9) where the minimum is taken over all flows 6 on Rf ut such that 

e 

div 9(x) = Afi 

for all x £ Ri. 


ii. £* = min £(60 where the minimum is taken over all flows 6 on R] n such that 

l e 1 

div0(x) = 4 fi(x) 

for all x £ Ri\di n Ri. 

With the preceeding preparation, the proof of Theorem 4.1 follows from a series of lemmas which we 
prove next. We first prove a monotonicity result. It is standard but we include it for completeness. 

Lemma 6.7. If a £ 11 satisfies Assumption 1 then for i = 1,2, 

£* < £ 
t'l < --i j 

where £i,£f are defined in (6.36) and (6.37) respectively. 

Proof. Wlog we take i = 1. By Lemma 6.3 and (6.30) we have £\ = Gi,n(x)Af\(x), and 

similarly by (6.33) we have £± = ^ x eR in n( x )^/i( x )- The lemma now follows since 

G*,n < Gi, n (6.39) 

which immediately follows from definitions of Gi >n , G\ n ((6.24) and (6.32)) and the simple fact that, 

riRiUdinlh) <r(Ri). 

□ 

The previous lemma shows that £i > £*. The next lemma says that nevertheless they are close to 
each other. 
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Lemma 6.8. There exists a f3 > 0 such that if a € Q satisfies Assumption 1 then for i = 1,2, 

where the constant in the O depends on 5 appearing in Setup 1. 

To prove the above we need a few preliminary lemmas first. Let, 

o* = V {Gin). (6.40) 

Thus (6.37) is the same as £( := £(#*). Hence by Remark 6.5, ii. has minimum energy among all 
the flows on Rf 1 with the same divergence condition. Even though Rf 1 is obtained from a subgraph 
of HJ^ by identifying some vertices we will denote the edges of this graph by the corresponding edges 
of U n . The next lemma is a conservation of flow result. It says that the total flow across di n R\ is 
at most a constant. 

Lemma 6.9. 

Y Y 

czf) R y ~ X 

XbOinit! 

Proof. We start with the observation that for all such x , y as in the sum, G\ n (x) = 0 and G\ n (y) > 
0 which follows by definition. Hence 6*{x,y) > 0 for all such ( x,y ). Now by (6.12) 

Yj div(0*)(z) = 0. 

z&R[ n 

By (6.33) for any z G R\ n such that z G R\\di n Ri 

div(0*)(z) = 4/i. 

Note above that even though the graph is not regular the 4 appears since every vertex in the 
support of f\ has degree 4. Also 

Y div(0*)(z) = Y 4 /i + Y H 9*(v,x)- 

z£R\ n zeR 1 \d in R 1 x&di n Ri (J R\ ye .Rp^a in h x 

The LHS is 0 and the first sum on the RHS is at most 4. Thus we are done. □ 

We state two more lemmas. The first result is a Beurling type estimate which says that for any 
connected subset A of with large enough diameter which is at a certain distance away from Ui in 
the probability that random walk started from a neighboring site of A, hits Ui iTl before hitting A 
decays as a power law in n. The result is standard on the whole lattice. The proof in our case with 
the necessary adaptions is provided in [16]. 

Lemma 6.10. [L6, Lemma 5.1] Fix c > 0. Consider A C 4J n be *—connected (Definition 5). Also 
suppose that mm(diam(A), d(Ui, A)) > c (here the diameter is in terms of euclidean distance not 
graph distance). Then for large n, for all such A, 

supP x (r(Ui in ) < t(A)) < 
x~A nP 

for some positive (3, C depending only on c and HJ. Here x ~ A means that x (f A and there exists 
y E A such that x is a neighbor of y. 

The next result is a basic isoperimetry inequality saying that the boundary of the set R\ is not 
small. 
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Lemma 6.11. There exists a c = c(S, U) such that 

mm(diam(di n Ri), diam(d ou tRi)) > c 

where 5 appears in Setup 1. 

Proof. Clearly it suffices to just show that 

diam(di n Ri) > c 

since every vertex in di n R\ has a neighbor in d ou tR\. Recall that by definition Ui )Tl C R\ and 
HJ 2 n C Ri . We first prove the lemma for the disc D. It is easy to notice that for the disc there 
exists a c = c(5) such that for any ball B of radius c, ©i and D 2 are connected by a lattice path in 
D n \B. Now since di n R\ separates Bi in ,B2,n this shows that diam(di n R\) > c. The general proof 
now follows by using the bi-Lipschitz nature of the conformal map cj) in (1.4). □ 

We are now ready to prove Lemma 6.8. Roughly the following will be our strategy: wlog we take 
i = 1. Recall 6* (6.40) is a flow on R™. We will construct a flow 6\ on from 0* by defining 
the value of the flow on the edges from di n R\ to d ou tRi (the extra edges in R° ut not in R™) and 
keeping the value of the flow on all other edges same as 9 *. 9\ will have the same divergence as 
VG i)Tl . Thus 

<£i< £{0{), (6.41) 

where the first inequality is by Lemma 6.7 and the second inequality by Remark 6.5 i. Then we 
would show that 

m) - n = o(T) 

for some (3 and thus complete the proof. Formally we do the following. 

Proof of Lemma 6.8. We construct a flow 9\ on R° ut : For each y in R\ C\d rn R.\ choose x <E d ou tR\ 
such that y ~ x (such an x exists by definition of di n R\). Let 

9i(y,x) = ^2 for eac]l C d in Ri 

zeRi 

z~y 

9i(y',x') = 0 for all other edges with y’ in d in R \, x e d ou tRi 
= 0* everywhere else. 

We first claim that 

£1 < £(9i). (6.42) 

By Remark 6.5 i. this will follow if we can show for all y £ R\ 

div(0i)(y) = d y (AG 1} n)(y) = 4/i (y), (6.43) 

where d y is the degree of y in R° ut . We know by definition that 

div(0*)(y) = d y (AG* hn )(y) = 4/i(y) for all y € Ri\d in Ri. 

Also by construction 0i = 0* on all the edges except the boundary edges of di n R\. Now by (6.8) 

d(U ll „,a i „fli)>(J 1 /4. (6.44) 

Thus f\ | a ^ = 0. Hence to verify (6.43) it suffices to show that 

div(0i) la in ^ nRl = 0. 
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Let y G di n R\ n R,\ . By construction we know that there exists exactly one x G d ou tR\ such that 

01 (y,x) 7^ 0. 

Also _ 

0i (y,x) = E 9*{z,y) = E 

z£Ri z£Ri 

zr^y z~y 

where the first equality is by definition and for the second inequality we use the fact that 9\ = 0* 
for all the edges in the sum by construction. Thus 


div(0i)(y) = E 0i (w,y) = 0, 


wER° ut 

w^y 


and (6.42) is verified. By (6.41) the proof of the lemma will be complete once we show that 


Now we claim that 


£{ 0 i)-e\ *<o(-g). 

n p 

sup „ *52 0 *(Vi z ) = 0 (~p)- 

y£d in Ri zefli 


(6.45) 


First notice that di n R\ satisfy the hypotheses of Lemma 6.10: 

. By (6.44), d(3 in i?i,Ui) ><Si/4. 

• The connectedness hypothesis is satisfied by Remark 6.2. 

• The diameter lower bound follows from Lemma 6.11. 
Thus by Lemma 6.10, for any y G di n R\ and z ~ y, 

1 


P*(r(Ui, n ) < r(d in R i)) = 0(^g 

where the constant in the O term depends on 5\ and ILJ. Since /i is positive only on Ui, n and 0 
everywhere else by (6.32) it follows that 


G* ln (z) < (r(Ui )Tl ) < T{d in R{)) sup G\ (x) 

v ' xeUi, n 


< (r(Ui )n ) < r(d in Ri)) sup G i, n (x). 

' X£V ltn 

The last inequality follows since for all x G U n , by (6.39), G\ n (x) < Gi >n (x). Thus by Remark 6.3 
we get, 

6*(v,z) = G* ltn {z)<0{±i), 

where the first equality follows from the fact that 0* = VG* n and G* ln (y) = 0. Since 0\ = 0 * on 
all but the boundary edges, 

m)-£(0*) = o( E (Ei 0 *^0)D 2 ) 

y£d in Ri z ~y 

( su P= Ei 6 '*( z ’0)D( E Ei^’0)D 

y&di^Rx z ~ y y£d in Ri z ~y 

0 


< 


< 
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where the last inequality is due to Lemma 6.9 and (6.45). Hence 

Si-st <£(e 1 )-£{e*) = o{\) 

nP 

where the first inequality is by (6.42). □ 

Remark 6.6. From Lemma 6.8 it follows that 

S-S 1 -S 2 >S-St-S 2 -0{\). (6.46) 

n p 

By (6.38) and the above remark, to prove Theorem 4.1, it suffices to prove a lower bound for the 
quantity 

£-£{-£*. 

We start with the following weaker version which shows it is always non negative. 


6.1.3. Weak version of Theorem f.l. 

Lemma 6.12. If a E Cl satisfies Assumption 1 then 

£ - £{ - £ 2 > 0 

Proof. Let us take the optimal flow 0 for £ on U n i.e. 

0 = VG n . 

For i = 1,2, let, 


(6.47) 


0i ■= 0 

6 * := 9 


RT 


RV 


and 


be the restrictions on R° ut and Rf 1 respectively. Now by (6.17) and Remark 6.4 it follows that, 

div(9i)(x ) = 4 fi for all x € Ri (6.48) 

div(6*)(x ) = 4 fi for all x £ Ri\di n Ri (6.49) 

for i = 1,2 where fi was defined in (6.14) and (6.15). Thus the flows 6( and 9* on R° ut and Rf 1 
have the same divergence as and VG* n respectively. Hence by Remark 6.5, 

Si <£&),£*<£$)■ 

Recall the definitions of from (6.8). Since R\ = U n \C'i j / c is connected and {R 2 U ^ 2 ,n} Q 
by definition R 2 C C\ j.. Thus the edge set of the graphs R\ n . Rf' 1 are disjoint. Therefore 

£ = £(9) > £(9 |^„) + £(9 = £(9l) + £(6 2 ) > £{ + £ 2 ■ (6.50) 

Hence we are done. □ 


However as already mentioned before, proof of Theorem 4.1 demands a stronger quantitative version 

wired 

of Lemma 6.12. We now proceed to obtain that. Let U n be the same graph as U n but with the 
following set of vertices 

di n R 1 l^J{cc £ R\ : a(x ) = 2} £ R\ : a(x) = 1} (6.51) 

glued. The points that are glued are all the red vertices in R± (which cause the blue random walk 

o wired 

to stop) and similarly all the blue vertices in U n \Ri. Now in U n the glued vertices act as a single 
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vertex. Let us denote it by tn. 

In the proof of Lemma 6.12 we took the flow VG n on U n and restricted it to Rf 1 and R? 2 ut ■ To 

wired 

prove Theorem 4.1 instead of restricting VG n . we will construct another flow on U n and then 

wired 

restrict it. We first construct the flow on U n . Recall / from (6.16). 

wired wired 

Lemma 6.13. There exists a unique flow 6 on U n with the following properties: 

i. 

wired 

div( 0 ) Iji^uRalVn = 4/ 

wired 

div( 9 )(•) = 0 for all other vertices except tu. 

ii. 

wired 

9 = min £(g) 

9 

wired wired 

where the infimum is taken over all flows g on U n which has the same divergence as 0 . 

wired 

Note that using the above data and (6.12), div( 6 )(ro) is determined exactly. 


wired 


Proof. For brevity let us call the set of flows g on U n satisfying the above divergence conditions 
as T. That the set is non empty follows since clearly by (6.17) the flow V(G* n ) \ wired is in this set. 


Now let 


U„ 


wired 

£ = inf £(g). 


(6.52) 


By standard arguments using compactness, see proof of [26, Theorem 9.10], there exists a unique 

wired 

minimizer 6 of the energy minimization problem. □ 

wired 

Remark 6.7. Let 0 be as in Lemma 6.13. The restrictions, 


wired 

wired 


9 i 

= 9 

1 R\ n 

wired 

wired 


e 2 

= 9 

\ftOUt 


under Assumption 1 satisfy the following divergence conditions 


wired 


div( 6 i)(x) = 4/i(x) for all x € Ri\di n Ri 

wired 

div( 0 2 )(y) = 4/ 2 (y) for all y € R 2 . 

Proof. Note that no vertex in either R\\di n Ri or i? 2 was in the set of glued vertices in (6.51). Now 

wired 

by Remark 6.4 for i = 1, 2, / 1^= f%. The remark thus follows since 9 by definition satisfies the 
divergence conditions. □ 

wired wired 

Hence by Remark 6.5 we have £\ < £{ 9 i), and £ 2 < £( 9 2 ). 
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Remark 6 . 8 . Similar to (6.50) we have 


£ = £( 9 )>£( 9 |#„) + £( 9 \ A o Ut ) = £( 9 l ) + £( 9 2 ) > £{ + £ 2 - 
Therefore we have 


i wired 

£-£i-£ 2 >£-£i -£ 2 -0(^) = £- £ + 


n 


wired 

£ -£*-£2 


- 0 (\' 

K nP' 


wired I 

> £~ £ -°U) 


n 


where the first inequality follows from (6.46). Hence to prove Theorem 4.1 by (6.38) it suffices to 
show that 

wired 

£ - £ >c (6.53) 

for some constant c = c(e). 


wired 

Notice that £ — £ is nothing but the drop in voltage when points at various voltages are glued 

wired 

together to get 4J n from U n . This is the well known Rayleigh’s monotonicity principle. For more 
details see [27, Chap 2], However for our purposes we need a quantitative version. To estimate 
such voltage drops we state and prove two technical lemmas. Let Un be a graph obtained from U n 
by gluing certain pairs of points. Let A be a set of k pairs of points ( 27 , z[),... (zj-, z' k ) which are 
among the pairs glued to obtain U n . Recall from (6.47), 9 = VG n . Consider the restricted flow 
V G n |u , (the flow that VG n induces on U n ). Let 9 denote the flow on U n such that 

m = inf {£(<?)}, (6.54) 

9 

where the infimum is taken over all flows g on U n which has the same divergence as 9 . Existence 

of 9 follows by a standard compactness argument as in the proof of [26, Theorem 9.10], (this has 
already been used in the proof of Lemma 6.13). 


Lemma 6.14. Let A C U n , 9 and 9 be as defined above. Suppose for each 1 < i < k, there exists 
simple paths in U n joining ( Zi , z') of length di then, 


£(9) < £(9) - 


Ei=l Gn{Zi) Gn(z[)Y 


D 2 Ell di 

where D is the maximum number of paths that intersect some edge. 


Thus the lemma measures the amount of voltage drop when points at different potentials are glued 
together. 


Proof. W.l.o.g. assume for all i, that G n (zi ) > G n (z! i ). Take a path joining z\ to z% of length di. 
Think of it as oriented towards z. t . Since on U n , z, L and z\ are glued, the path becomes a directed 
cycle. Let us denote it by 7 j. Also let 7 j be the same cycle in the reverse orientation. We first 
create a new flow 9 a on U n by sending an additional amount of flow fi along 7 j for alH = 1 ,... ,k. 
Thus 

k k 

9 A (e) = 9(e) + ^ /31(e G 7 i) + ^ -/31(e G 7 )) 
i= 1 j =1 
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for all directed edges e, where every undirected edge in HJ n appears with both orientations. Now 
the additional flow along an edge is at most Df3 since by hypothesis every edge is in at most D 
many cycles. An easy computation by expanding the squares now shows that 

k k 

£(0a) < £{B) + 2/3 ^2(G n ( Zi ) - G n (z')) + D 2 0 2 ^ d*. 

2—1 2 

Optimizing the RHS over (3 we see that, 

inf C (9A ^ c (f )'I E*= l( G n( z i) - Gn(z'i ))] 2 

mf £( ^) < m -^7-• 

Notice that since all the paths 7 ,; become cycles in U n and sending additional mass along cycles 
does not change the divergence of the flow, div (@ y 4 |fj ) = div($ ||j ), (note however that 6a and 0 
do not have the same divergence on HJ n ). Hence by (6.54), 8(0) < £(9a lu ) < £(9a) and the proof 
is complete. □ 

By the above discussion and (6.53), to use Lemma 6.14 in the proof of Theorem 4.1, one has to 
estimate the df s and D appearing in the statement of the lemma. That is the purpose of the next 
lemma. The statement of the lemma is rather crude but will suffice for our purposes. 

Recall the definition of LJ^ from (4.5). 

Lemma 6.15. Given a smooth domain HJ for any small enough d > 0 (depending on V), there 
exists D,a > 0 such that for all 5 < 5o(d) and all n = 2 m > N(5) the following is true: For 
any set A C U n containing dn points in distinct IjW ’s , there exists at least dn /10 disjoint pairs 
(zi, z [),... (z dn / 10 , z' dn/ 10 ) of elements of A such that: 

i. \G n (zj ) - G n (z'j)\ > a for j = 1,..., dn/10, 

ii. There exist paths in HJ n between Zj and z'j of length at most Dn, such that no edge in HJ^ 
intersects more than 2 paths. 

Remark 6.9. It will be useful later, to note that any set A C HJ n with |A| > cn 2 for some constant 
c independent of n satisfies the hypothesis of the above lemma for a certain value of d = d(c, U). 

Proof of Lemma 6.15 ii. Recall from (4.5) that all the U”’s have width at least 100/n. Now let 
Zi G dn . We will pair Zi from different ’s Now the only difficulty can come if the points he 
close to the boundary where the graph does not look like Z 2 . However in that case for each such 
Zi one can take a path in HJ n away from the boundary of the set HJ, ending at say yi, which lies 
entirely in U^ -1 ^ U U U^ +1 \ (smoothness of LJ (see (9.1)) allows us to do that. Now to avoid 
intersection of these paths we will only consider points z* in HJ® that are not consecutive. 

The proof is now completed by pairing points y* into pairs ( Vj,y'j) such that corresponding ILJW for 
the elements of a pair are at least d/10 apart. That the elements of each such pair ( yj,y '•) can be 
joined using disjoint lattice paths of length at most Dn for some D depending only on U is easy to 
see. We omit the details. See Fig 10. 

Proof of i. Recall that by the above construction for any pair Zj and z'j the corresponding HJ^’s are 
at least d/10 apart. Thus by Theorem 5.1 and Lemma 5.3 there exists do such that given 6 < do 

lim inf [inf (G n (zj ) - G n (z'A\ > c, 

n=2 m j J 

m—>00 

for some constant c depending only on d and U. □ 
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Figure 10. Figure illustrating the proof of Lemma 6.15 for the case of the circle 
where points in different U® are connected by almost edge disjoint paths. 

6.2. Proof of Theorem 4.1. In combination with Theorem 6.1 it suffices to prove Theorem 4.1 
under Assumption 1. To this end we will show: for all a E there is a set A which is a subset 

wired 

of the set in (6.51) (hence is glued to get U n from U n ) such that: A satisfies the hypothesis of 
Lemma 6.15. 

Before showing the above we first discuss why it suffices. Immediately using Lemma 6.15 we have 
pairs (zi, z[),... (zk, z' k ) with k = H(n) which are glued and have paths connecting them of length 
0(n ) such that no edge appears in more than 2 paths. Also, G n (zj) — G n (z'j ) > d , for j = 1... k. 

wired wired 

Recall £ from (6.52). Thus using Lemma 6.14 we get, £ < £ — d, for some constant d = d{e). 

Hence we are done by (6.53). 

We now proceed to showing existence of such sets A (the arguments will be rather crude and not op¬ 
timal in any sense. However since we are interested only in statements up to constants this suffices): 

Recall P from (1.6) Let 

P-4 > P-3 > P-2 > P-i > Po = P > P\ > Pi> P:i > P 4 (6.55) 

such that for all i = —4, —3,..., 3, 

d{iPnlPi+i) = e i ( £ i is specified soon). 

See Fig. 11. Recall (6.8) and (6.9). Now consider the following cases: 

i. dinRi is a subset of Ug_ 2 , 

ii. di n R\ is a subset of U^ 2 , 

in. di n R\ intersects U J g. i and HJ^ for some i E {—4,... ,3}, 
iv. di n Ri is a subset of U^_ 3 FI Ug 3 , 

We choose e± < ^ such that, 

area(U^._ 1 flU^) < e/100 (6.56) 

for all j = —4,... 3. (100 is just a large enough number and is nothing special) Also let e' > 0 be 
such that the left hand side in (6.56) is at least e'. 

The arguments for cases i., ii. are symmetric and hence to avoid repetition we provide arguments 
only for case i. 
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Figure 11. 73 .’s on either side of 7 ^ to measure how far the boundary of R \, the 
blue curve , deviates from 7 ^. 

By hypothesis di n R\ lies entirely in Ug_ 2 . Now by arguments similar to the proof of Lemma 6.15 
ii. for large enough n, any two points in U n nU^ j are connected by a path in U n which lies entirely 

0 0 C 

in U n nU^ i+ei / 2 (hence does not intersect di n R\). Thus U n flU^ is in R\ . Recall by definition 
area(U^ o ) = (1 — a)area(U). Thus, area(U^_ i ) > (1 — a)area(U) + e' . Hence number of blue vertices 
in U n n is atleast e'n 2 . Recall from (6.51) that both (j _ 1 ( 2 ) n R\ and (j _ 1 (l) n R\ are glued 

wired 

to obtain HJ n from U n . Thus at least e'n 2 vertices are glued and hence we are done by Remark 6.9. 

Case in. By hypothesis in this case di n R\ (which is a *— connected set by Remark 6.2) intersects 
different level sets of the function G n . Thus the intersection with the different level sets could be 
taken to be A in the usage of Lemma 6.15. 

Formally we do the following: by hypothesis di n R\ intersects both Ua ; and UJg. | Now by compact¬ 
ness there is some constant c = c(/3, ei, U) such that for all i = —4,..., 3, 

d{ 7/3 i ,7ft+i) > c 

Since d in R\ is *—connected (Remark 6.2), it intersects at least dn many U^’s for some d = d(c). 
Thus A = di n Ri satisfies the hypothesis of Lemma 6.15 and is also glued by (6.51). Therefore we 
are done. 

Case iv. Recall that by hypothesis a G IF-j and di n R\ is a subset of U ^_ 3 nUg 3 . Since a £ QC and 
e i < jog! there exists at least dn many HJ^’s containing at least one red vertex all of which are in 
U n nU ) g _3 or similarly there exists at least dn many HJ^’s containing at least one blue vertex all of 
which are in U n H U| . We take the set of such points to be A. 

The proof of Theorem 4.1 is thus complete. □ 


7. Proof of Theorem 1.1 

Given the preceding results the proof of the main theorem follows by similar arguments as in the 
proof of the main result in [ ]. We first state the following two hitting time results. Recall the 

sets from Table 1. 
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Lemma 7.1. Given any e > 0 there exist a positive constant e\ such that for all small enough 
5 (blob size) there exist positive constants c = c(e,S),d = d(e,5), such that for all large enough 
n = 2 m and a 6 n A ei 

P CT (r(<? e i/4) > dn 2 ) < e~ cnl/3 . (7.1) 

Lemma 7.2. Given e > 0 there exist a positive constant e' such that for all small enough 6 (blob 
size), there exists b = 6(e, S),d = d(e, 5) > 0 such that for all large enough n = 2 m and a e G e FF e / / 2 

K(T(g c el/4 ) > e dnl/3 ) > 1 - e-^ 3 . 

The proof of the above two lemmas involve arguments involving IDLA on HJ n and is postponed to 
the appendix. We now show how to finish the proof of Theorem 1.1 using the above. 

Lemma 7.3. Let r e = inf{f : S(t) G Q e }. For any e > 0, for small enough 5, there exist positive 
constants c = c(e, 5),d = d(e, 6), N = N(6) such that for all n = 2 m > N and all subsets S C U n of 
cardinality [an 2 J we have 

P(r £ > dn 2 | 5(0) = S) < e~ cnl/3 . 

Proof. We will prove something stronger which will be used in the proof of Theorem 1.1. However 
for brevity we first need some notation. We start by recalling the sets in Table 1. Given any 
positive 6\ and e 2 define the set 

Cei,e2 = Gei G r e2 . 

We claim that for small enough e, e' > 0 there exists constants D,D' > 0 such that for all large 
enough n = 2 m and any cr G Ft, 

IV('r(C £ VV/2 )>Dn 2 )<e- D ' nl/3 (7.2) 

E 4^(C e i/V/ 2 ) ^ 2Dn 2 ) (7.3) 

(7.3) clearly follows from (7.2) as the latter holds for any starting a and hence the expected hitting 
time is clearly dominated by Dn 2 times a geometric variable with mean 2. Also showing (7.2) 
clearly proves Theorem 7.3 since C e i/4 e y 2 C Q t 1 / 4 . Let 

t' = r(T e , /4 ),T" = t(T‘ /2 ) 
t'" = t(H( £) n A £i ),t"" = t(Q £ 1/4). 

where the ei in the definition of t'" is the same as the one in the hypothesis of Lemma 7.1. Now 
to show (7.2) we notice the following containment of events for any positive A: 

[•{V < An 2 } n {t" > e cn | a T f] n {t"' < An 2 \ a T '} n {t"" < An 2 \ <r T w}] C {v(C e i/4 )e // 2 ) < 3 An 2 }. 

To see why this containment holds we first notice that T e / / 4 C T e // 2 . Thus the first two events imply 
that the process hits the set T e // 2 in An 2 steps and stays inside for an exponential (in n) amount of 
time, and in particular stays in T e // 2 from time An 2 through time 3 An 2 . In addition, the third and 
fourth events together imply that regardless of where the chain is at time An 2 , the chain enters 
n A ei by time 2An 2 and then enters f/ e i/4 by time 3 An 2 . Hence in the intersection of the four 
events, the hitting time of C e i/4 e ,/ 2 = ^ e i/ 4 H T 6 // 2 is at most 3An 2 . Let D = 3 A. Now for a large 
enough constant A, there exists D' > 0 such that the probabilities of of all the four events on the 
left hand side are at least 1 — e~ D ' nl/3 . This follows by Lemmas 4.4, 4.5, and 7.1 respectively. 
Hence by the union bound, (7.2) follows for a slightly smaller value of D'. □ 
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7.1. Proof of Theorem 1.1. Note that (7.2) is true for all small enough choices of e and e'. 
However given e > 0, by Lemma 7.2 there exist e',b,d > 0 such that for all large enough n = 2 m , 
and a 6 C e e i / 2 

IMt(<? £ c 1/4 ) > e dnl/3 ) > 1 - e~ bnl/3 . (7.4) 

For such an e' the proof of Theorem 1.1 follows by using Lemma 4.7 with the following choices of 
parameters: 

A = C e , e7 2 , B = g c eX/ 4 , H = 2Dn 2 , t 2 = e dnl/3 

The above choice of parameters satisfy the hypotheses of Lemma 4.7 by (7.2) and (7.4). □ 

The remainder of the article provides a sketch of the proof of Theorem 5.1 and finally in the 
appendix we finish off with the proofs of Lemmas 7.1 and 7.2. 


8. Outline of the Proof of Theorem 5.1 

We provide the sketch of the main steps in the proof of Theorem 5.1. The full proof appears in 
[16]. The proof relies on results based on convergence of random walk on Un, to Reflected Brownian 
motion on HJ. For a formal definition of Reflected Brownian motion on U see [11, Definition 2.7], 
[5, 10, 9]. Throughout the rest of the section we will denote it by Bt. The basic idea behind proving 
Theorem 5.1 is to show that the function G*(z) is the same as, 


G(z) = 


area(Ui 


[P z (B t E Ur) — P z {B t E U 2 )] dt — c. 


( 8 . 1 ) 


The constant c is chosen such that the integral of G along <9HJ is 0. Formally, we parametrize the 
boundary dU by 0 E [0, 27t) via the conformal map 4> (1-4). We then fix c such that 


[ Go <XC)^ = o. 
'ICI=1 C 


( 8 . 2 ) 


Compare the above expression with G n (-) defined in (4.1) (recall that the constant c is the same 
in both expressions). The proof of Theorem 5.1 has two parts. One part shows that G n converges 
to G, (see Lemma 8.2). This will follow by convergence of the random walk measure on HJ n to Bt. 
For more on this see [9], [ I ] and the references therein. The remaining step then is to show that 
indeed G* = G. This will be proved using the fact that the density for Reflected Brownian motion 
is a fundamental solution to the Neumann problem and hence the function G is roughly satisfies, 


d 2 d 2 \ _ 4 

dx 2 dy 2 ) area(Ui) 


1(U 2 ) - l(Ur)) 


(8.3) 


with Neumann boundary conditions. One then checks that G* is a solution to the above problem. 
The proof is then complete by uniqueness of the solution of such a problem. For brevity we will 
adopt the following notation 

d 2 d 2 


A := 


dx 2 dy 2 


Recall (5.3). Thus we use A to denote the laplacian in both the continuous and discrete setting 
since there will be no scope of confusion. 
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8.1. Closed form of the limit. 

Theorem 8.1. [16, Theorem 4.2] Let U be as in Setup 1. For all z G U, 

G(z) = G*(z)- (8.4) 

As discussed above, to show this informally one identifies G as a solution to a second order partial 
differential equation also satisfied by G*. The result then follows by uniqueness of such a solution. 
To give a general idea we quote a result in the theory of boundary value problems with Neumann 
boundary condition which is used in the proof. 

Lemma 8.1. Let function w G C' 2 (U)nC' 1 (U) be a function on U satisfying the following properties 

A (w) 

dw 
dv 


f 

4 

= 0 on <9U 


where denotes the normal derivative and f G L 1 (M,U). Then 

wo <f>(z) = d+- [ f o 0(C) I ^ 1 2 (C) log(|(C — z)(l — (z)\ 2 )dfdr], (8.5) 

where 

j-f ^m d(=d . 

2 TTl Jq 0 C 

Proof. The result for the case U = D (the unit disc) is stated as [6, Theorem 8]. The above lemma 
now follows by change of variable under composing with the map <j>. □ 


Theorem 8.1 now follows roughly by arguing that the function G(-) satisfies Lemma 8.1 with 
AG = area m 1 ) (1(U 2 ) — 1(U2)) and noticing that G* is nothing but / in Lemma 8.1 replaced by 

areafUi) — 1(^2))- However formally one works with “smoother” versions of (1(U2) — 1(U2)). 

The detailed proof appears in [16]. 


Remark 8.1. Recall the constant c in (8.1) such that. 

[ G o = 0. 

j ICI=1 C 

Consider the the special case of the unit disc i.e. U = D with point x\ = —i, X 2 = i and y\ = 
— (1 — 6 )i, y2 = (1 — S)i. Now owing to the invariance of Reflected Brownian motion on D under 
the transformation £ —> —one sees that the integral in (8.1) 

2 r°° 

-— / [P z (B t G Dr) - P z (B t G D 2 )] dt 

area(Di) J 0 

is odd, (because with the choice of yi,y 2 , we have Di = —D 2 ). Thus c is 0 and 

2 r°° 

G(z) = -— / [P z (B t G Di) - P z (B t G D 2 )] dt. 

area(Di) J 0 
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(xi,X 2 ) . 



Figure 12. i. Locally the region near the boundary looks like a half plane, ii. 
Illustrating the proof of Corollary 9.1. 


8.2. Convergence. The last missing piece is to show that the function G n converges to G. To this 
end we now state the following lemma which along with Theorem 8.1 proves Theorem 5.1. 

Lemma 8.2. [16, Lemma 4.8] For all small enough 5 , 

Jim sup | G n (x) - G(s)| = 0 

71=2™ a:eU 

where G is defined in (8.1). 

Note again the suppressing of 5 in the above statement. The proof of the lemma uses the local CLT 
estimates for random walk approximation of Reflected Brownian motion obtained in [11, Theorem 
2.12], The complete proof appears in [16]. 

Proof of Theorem 5.1. The proof follows immediately from Lemma 8.2 and Theorem 8.1. □ 

9. Basic geometric properties of U 

We include a short discussion of some basic geometric properties of U n used before in the article. 
The following is a well known property of the boundary of U as in Setup 1: since the boundary 
9U is analytic there exists aC>0 and an eo such that for all x G 9U there exists an orthogonal 
system of coordinates centered at x = (xi,x 2 ) such that for all e < eo 

B(x, e) n U = { (x\ , xj) G B(x, e) : x\ G (xq — e, x\ + e), x' 2 > /(aj)} (9.1) 

and 

| f{x[) - x 2 \ < C\x\ - xq| 2 . 

The above is a simple consequence of Taylor expansion up to second order of the curve locally near 
x. See Fig 12 i. As a simple corollary of the above fact we see that HJ satisfies the following property 
which shows that the yfis in Setup 1 can indeed be chosen. 

Corollary 9.1. Let HJ be as in Setup 1. Then there exists So = <5o(U) such that for all x G U and 
5 < 5 q there exists y G U such that d(y , x) < 5 and B(y,S/ 2) C U. 

Recall that d(y , x) is the euclidean distance between x and y. B(y, 6) denotes the euclidean ball of 
radius 6 with center at y. 
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Proof. Choose 60 < eo/4 such that C 6 q < where eo and C appear in (9.1). For any 5 < 60 the 
lemma is immediate if d(x,<9U) > 6 / 2 . since then we can choose y = x. 

Otherwise let z = (z\,Z 2 ) G dU be the closest point on the boundary to x. Now in the local 
coordinate system centered at z as in (9.1) choose y = (z\, Z 2 + 6/2). Then 

d(x, y) < d(x , z) + d(z, y) < 6 . 

Also clearly B(y , 6/2) C U and hence we are done. See Fig 12 ii. 

□ 


Concluding remarks 

Note that Theorem 1.1 only proves concentration of the interface. 

Fluctuations. A natural next step would be to find the order of magnitude of fluctuations of the 
interface. 

Higher dimensions. One can also study the process in higher dimensions and try to characterize 
the hyper surface separating the two colors at stationarity. It is natural to guess that it should 
be the level set of the analogous higher dimensional potential function which is harmonic in the 
interior of the corresponding domain and has positive and negative singularities at the two sources 
with Neumann boundary conditions. 

Multiple particles. Proving a theorem analogous to Theorem 1.1 for Competitive Erosion with 
more than two particles remains an intriguing challenge. 

Appendix 

We finish off with the proofs of Lemmas 7.1 and 7.2 that were used in Section 7. We start with 
a generalized IDLA estimate. For definitions and other results on standard IDLA on the square 
lattice see [ ]. However our setting will be different in two aspects: 

• The graph will be HJ n instead of Z 2 . 

• At time 0 some of the sites will be filled unlike the standard setting where every site is 
vacant initially. 

Formally we have the following: recall HJ^ from (4.5) and U( Q ) from (1.7). We start an IDLA 
process on U n with the random walks starting uniformly from U 2 (see Setup 1) with the following 
initial condition: All the sites in U n \{HJ n nU( Q )} are filled along with an arbitrary e'n many U^’s 
for some small e' > 0. Let us now fix an e > De' ( D is just a large enough constant to be specified 
later). The next result will bound the size of the cluster after en 2 particles have been released 
similar to standard bounds (see for eg: [ ]). Let the IDLA cluster (set of filled sites) at time t be 

A t . 

Theorem 9.2. (IDLA upper bound) There exists c, C depending only on HJ, a such that given small 
enough e > 0 and starting with any initial condition as described above, for all large enough n, with 
probability at least 1 — e~ c W en ) } 

A en 2\A 0 C [{U\U( a _Cv^)} n U n ]. 

Remark 9.1. Note that the distance between the geodesics which act as boundaries of U( a ) and 
U (Q _ e) are up to universal constants, e apart for all small e. Thus the above theorem is rather crude 
and says that the aggregate of en 2 particles cannot grow beyond an ‘annulus’ of width \fe (recall 
that the initial cluster was inside U\U( Q ) except for e'n many U^’s) even though one should expect 
it to be e. However it suffices for our purpose. 
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Figure 13. 71 and 72 denote the boundaries of U( Q ) and U ( a -Cy/e) respectively. If 
the intial aggregate contains all of the region above 71 i.e. (U\U( Q )) and at most 
e'n where e r <C e many thin strips (see (4.5)) then in en 2 rounds of IDLA on 
U n where the random walks start from U 2 ,n, the aggregate does not occupy any site 
below 72 i.e. the new part of the aggregate stays in {U\U( a _ c .^}. 

The precise definition of IDLA is provided along with the proof of the above. 

9.1. Proofs of Lemmas 7.1 and 7.2. With the help of the above theorem we now finish all 
the remaining proofs. A corresponding lower bound analogous to Theorem 9.2 will be the content 
of the proof of Lemma 7.1. We follow the approach in [ ]. Similar statements for IDLA on the 

entire lattice Z 2 have cleaner proofs owing to symmetry of Z 2 . However in our setting we need more 
robust approaches and we will rely on certain heat kernel bounds for the random walk on U n from 
[11], which are summarized in [16]. We first quote the following hitting time estimate. 

Lemma 9.1. [16, Lemma 5.3] Given small enough e > 0, for all z E U n n U/ a _^, and y E U n 
such that d(y, z) < e 2 , 

F y {r(z) < r(U n \{U n nU (a _ e) })) = 
where the constant in the 0 notation depend on e, a, ILJ. 

Proof of Lemma 7.1. Since by hypothesis a E fL e ) r\A ei it has at most ein 2 red vertices in U( Q ) nHJ„ 
and similarly at most ein 2 blue vertices in UW nD„. To show that the configuration reaches G e i /4 
we need to show that the red particles in U a _ e i/4 gets wiped off and similarly the blue particles on 
the other side. Because of the obvious similarity in the situations we will provide argument only 
for the former case. 

Now notice that the red cluster in erosion is clearly dominated by the IDLA process as described in 
Theorem 9.2 with the initial aggregate being the set of red vertices in a. Thus by Theorem 9.2 with 
probability at least 1 — e c (%/ fn ) no red walker reaches ^f a -c 2 ^ in the next C\en 2 rounds where 
C *2 depends on C\. Observe that on the above event the blue cluster insider U( a _c? 2V ^ dominates 
a ‘killed’ IDLA process where the blue walkers start uniformly from Ui )n and perform IDLA but 
are killed on hitting the boundary of the set U n n U r a _c 2y /e) (since there is no interaction with 
red random walkers). Thus it suffices to show that in the next Cien 2 rounds of the ‘killed’ IDLA 
process the blue walkers will occupy and hence wipe off all the red vertices in ILL a _ £ i/ 4 ) (at most 
ein 2 ), which amounts to proving a lower bound for the ‘killed’ IDLA aggregate for the blue walkers. 
We adapt the proof in [24]. For any positive integer i we associate the following stopping times to 
the i th blue walker: 

• a 1 : the stopping time in the killed IDLA process 
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• rl: the hitting time of z £ U n 

• t 1 : the exit time of the set U( Q . -Ciy/t) nU ™- 

For 2 6D„n U( a _ e i/4) we define the random variables, 

Ciera 2 

N = 1 {ri<a i )i the number of blue walkers that visit z before stopping 

2=1 
C\en 2 

M = l( r i <T i), the number of blue walkers that visit z before exiting U( a _c? 2v ^) 

2=1 

Cien 2 

L = l(ir i < T| < T i )i the number of blue walkers that visit z before exiting U ^ a _c 2y /i) 

2=1 

but after stopping. 

Note that random walkers can stop in two ways: either by occupying an ‘empty’ (red) vertex or by 
exiting U ( a -c 2 ^/e)- Thus N > M — L. Hence, 

P (z$ A c i en2 ) = P (N = 0) < P(M < a) + P(L > a). (9.2) 

where the last inequality holds for any a. Now by definition, 

E(M) = C ie n 2 P Ulin (r(z) < r(U n D {U\U (Q _ CaVi) })), (9.3) 

where Pui„( - ) denotes the measure when the random walk starts uniformly from Ui n . We now 
bound the expectation of L. Define the following quantity: let independent random walks start 
from each ‘empty’ (red) site in HJ r a -c 2 ,/e) an h let 

L = ^ l( r (- 2; ) < T(U n n {U\U( a _c< 2v ^}) for the walker starting at w). 

w ^(a-C 2 Ve) nU " 
w-.empty 

Clearly L is stochastically dominated by L. Hence the RHS of (9.2) can be upper bounded by 
P(M < a) + P(L > a). Now 

E(L) = ^ P w (t(z) < t(U n C {U\U (a _ C2Vi) })). 

«jeu (a _ c . 2v ^ ) nu„ 

w-.empty 

Note that by hypothesis there exists at most ein 2 empty (red) vertices in E^ a _ C2 ^ n U n . Since 
the above sum is over a set of size at most ein 2 , by Lemma 9.1 we have the following, 

E(L) < 

log n 
n 2 

E (M) > C"C\€- - 

logn 

for some constants C 1 , C". The first bound follows by summing the bound in Lemma 9.1 over a set 
of points y , of size ein 2 (this is maximized when all the points are in a ball of radius yTi around z). 
The latter bound follows since uniformly for any point z £ U n n U( a _ £ i/4), the blue random walk 
starting uniformly from Ui n has a chance Ci = C(e, U) of reaching a point e 2 close to z before 
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exiting U t a _c 2 ^/iy l Now once the random walker reaches such a point, using Lemma 9.1 we get the 
bound of C"i We now choose ei small enough so that E(L) < E(M)/4. Choose a = E(M)/2. 
Now using Azuma’s inequality for indicator variables we get 

P (L > a) < exp (—dn) and P(M < a) < exp(— dn) 
for some constant d > 0. Thus in (9.2) we get 

P(M < a) + P (L > a) < 2 exp {—dn). 

The proof of the lower bound now follows by taking the union bound of the at most ein 2 , empty 
2 G U n n U( a _ e i/ 4 ) i.e. such that a(z) is red. □ 

9.2. Proof of Lemma 7.2. For any set A C define the positive return time t + (A) to be 
inf{i > 1 : at G A}. We break the proof into a couple of lemmas. Recall Table 1. 

Lemma 9.2. Given positive numbers e,c there exist e' such that for all small enough 5 (blob size) 
there exists D = D(e, c, 5) such that for any a G n T e t 

P CT (T + (fi (e )) > cn 2 ) < e~ Dn2 , (9.4) 

for all large enough n = 2 m . 

Proof. The proof follows from Lemma 4.6 i. First let a' = (J\ for the competitive erosion chain 
starting from ao = a. If a' G fL e ) we are done. Otherwise since a G r e /, 

x := w(a') > w m ax - e'n 2 . 

Now let 

X{t) = w(a t+ 1 ) 

where w(-) is the weight function defined in (4.3). Thus X(0) = x. We make the following choice 
of parameters: 

B = fl( e ), Ai = w max , A 2 = log(l/<5) 
a\ = o(e) appearing in Theorem 4.1 
a 2 = e'n 2 with e' < a\c/2 , T = cn 2 . 

The choice of a 2 works since a G T e / by hypothesis. Also the choice of e' ensures that a 2 — a\T < 0 
and hence the hypothesis of Lemma 4.6 i. is satisfied. Thus for e' = Ljr by Lemma 4.6 i. for all 
a g nr e /, 

F’ ( 7 '(' r (^L)) >T) < exp (-Dn 2 ) , 

for some D = H(e,c, d). Since r + (0( e )) starting from a is one more than r(fl( e )) starting from a' 
we are done. □ 

Lemma 9.3. For all small e > 0 there exists e' > 0 such that for all small enough 5 (blob size) 
there exists a constant b(S, e) such that for large enough n = 2 m and any a G f2( e ) 0 r e / 0 Q e 1 / 4 , 

P CT (r+(L! (e) ) < r(fi\a el/4 )) > 1 - e ~ bnl/3 . 

1 Formally once can take a tube of a small enough width joining a point in the ‘interior of’ U n (says lOe away from 
the boundary) to any point 2 and see that the random walk stays in that tube and hits the end with the point z 
with constant probability Ci = Ci(e, U), independent of n. That the random walk starting from Ui,„ reaches such 

an interior point before exiting P( a -c 2y fc) is a straightforward consequence of the heat kernel estimates [16, Theorem 
5.2], 



COMPETITIVE EROSION IS CONFORMALLY INVARIANT 


49 


Proof. Since a £ n Q t i/ 4 , by Theorem 9.2 (since clearly IDLA dominates erosion), 

P CT (T(ft\g e i /4 ) > dn 2 ) > 1 - e~ hnl/3 

for some d, h depending on e, <5. Now using c = d/2 in Lemma 9.2 we get that we can choose e 7 such 
that for any o £ fb e ) n T e / n C? e i/ 4 

F a (T + {n (e) ) > dn 2 / 2) < e~ Dn \ (9.5) 

for some positive constant D. Thus for such an e 7 , 

P ff (r+(n (e) ) < r(n\g £ i /4 )) > P CT (r + (fi (e) ) < dn 2 / 2) - P„(T(fi\G e i /4 ) < dn 2 ) 

> l - e - hnl/3 - e~ Dn2 . 

Hence we are done by choosing b = □ 

9.3. Proof of Lemma 7.2. We will specify e 7 later. However notice that for any small enough e 7 
if <7 £ r e72 then by Lemma 4.5 there exist positive constants c, d depending on e 7 and the blob size 
such that for large enough n = 2 m , 

P cr (r 7 > e cn2 ) > 1 - e~ dn2 (9.6) 

where r 7 = r(fi\r e /). Now notice that by hypothesis a £ Q e and hence a £ fi( ce ), for some constant 
c = c(U). Let t ( 1 ), t( 2 ) ... be successive return times to i.e. t 7) = 0 and for all i > 0, 

T(i+1) = inf {^ : t > T W , a t £ 0 (ce) }. 

The following containment is true for any positive b Let s := . Then 

{r(g/ i/ 4 ) < s} C < s such that a t (f T e /}U | [J |{r 7) < t{Q/ 1/4 ) < r 7 + i)} n {a T(j) £ T e /}| 

V'=° 

The above follows by hrst observing whether there exists any i < e b ' nl/3 such that < 7 j ^ T e /. Also 
let j be the hrst index such that 

T U) < T (^ e C i/4) < r (i+i)- 

Notice that on the event {t(Q c 1/4 ) < s}, we have j < s. Also by definition on the event {cr* £ 
T e / for all i < sj, we have a T(/) £ T e / n fl( ce ) n C? e i/ 4 . Thus by the union bound, 

P a (r(g/ 1/4 ) <e bn/ ) < P(3i < s such that < 7 j ^ T e /) 

e b'n 1/3 

+ Y Pct ({°>oo e r e / n n {ce) n g e i /4 } n {r(g e c 1/4 ) < t u+1) }). 

3 = 1 

Recall that by hypothesis a £ T e // 2 - Thus by (9.6) for any small enough e 7 there exists d = d(e', 6) > 
0 (where S is the blob size) such that for any b' > 0 the hrst term is less than e~ dn2 for large enough 
n. Also notice that by Lemma 9.3 we can choose e 7 such that every term in the sum is at most 
e -bn 1/i f or some constant b > 0. Thus for such an e 7 , putting everything together we get that for 
any b' > 0 and large enough n, 

ei /„ 1/3 

P CT (r(S e C i/ 4 ) < e b ’ nl/Z ) < e~ dn2 + Y ■ 

Z=1 
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Hence by choosing b' < b we get that for large enough n, 

P(t(£ £ c 1/4 ) < e b ' nl/3 ) < 2e( b '~V nl/3 . 

The proof is thus complete. □ 

All that is left to be done is to provide the proof of Theorem 9.2. The argument will be a modifi¬ 
cation of the argument in [12] adapted to our setting (for the sake of reference we will follow the 
same notation as in [12] throughout this section). Let S C U n . We now precisely define the IDLA 
process in our setting. In order to define IDLA, we first define adding one particle started at x to 
an existing aggregate S. For x G U n , denote by A(S\x) the IDLA aggregate obtained as follows: 
Let £ = (£(0),£(1),...) be a random walk on Un started at £(0) = x and let ts be the first time 
this walk is not in S. Define A(S\x) := S U {£(ts)}- Also consider a slightly more general process, 
where the growth of the aggregate is stopped at certain stopping times, e.g., upon exiting a set T. 
Denote by A(S;x —» T), the aggregate obtained by letting a particle randomly walk from x, but 
pausing it if it exits T, which is defined as follows. Let £ be a random walk on U n started at x. Let 
ts be the first time this walk is not in S as above, and let tx be the first time £ exits T. Define 
A(S; x —» T) := S U {£(ts A (tx — 1))}- To keep track of the position of a paused particle, define 

p( g. x > rp\ _ f £(Lr) if tx < ts, 

' ’ { _L otherwise, 

so _L means that the particle is already absorbed in the aggregate. Given vertices x±,...,Xk in 
U n and a set T, define A(S\x\, ...,Xk —> T) to be the IDLA aggregate formed from an existing 
aggregate S by k particles, started at x\,X 2 , ■ ■ ■ Xfc and paused upon exiting T. Formally we use 
induction: Sq = S, Sj = A(Sj-i) Xj -» T ) for j G {1,2,... A:} and A(S] x ±,..., Xk —> T) = Sk- 
Again to keep track of paused particles, define P(S; x \,..., Xk —> T) to be the sequence of particles 
paused in this process, i.e. if pj = P(Sj- 1 ; Xj -A T) for j G {1,..., k}, then P(S‘ x±, ..., Xk —> T) is 
the sequence (pj : pj ^_L). We recall the Abelian property of IDLA: A(S; xi ,..., Xk) has the same 
distribution as A(A(S; x\.,Xk T); P(S‘ x ±,..., Xk —> T)). For more details about the Abelian 
property see [24] and the references therein. 

Proof of Theorem 9.2. We begin by stating two lemmas from [ ] adapted to our setting. Recall 

from (1.5), (should not be confused with U n which denotes the graph) and the ‘shells’ from 
(4.5). 

Lemma 9.4. [12, Lemma 5] Fix an interval IgM. There exists small constants e, e± such that for 
all large enough n and positive integer r with n 1//6 < r < en and /3 G I, the following holds: Let 
x G U n PI {U\Ug} and let S C {U\Ua + i:} H U n be so that [5 PI [U n PI Ug]| < eir 2 . Moreover let A 

be a set which intersects at most e± fraction of the U^ ’s between U^ and U^+n, Let £ be a random 
walk started at x and stopped upon hitting Un + 1 . Then, 

P(£n (u 0 \(S u A u v p+ r _)) + 0 ) > v , 

for p = p(I, U). 

Recall U^/ C U^ for any f3' > /3. The above lemma says that if the size of the intersection of the 
set S with the region between U^ and U^_|_r is not large then there is a constant chance that the 
random walk starting from U n n{U\Ug} will visit sites in U / g\Ufl + r which are neither in S nor in A 
before hitting Ug . r. After analyzing the behavior of a single particle, we can analyze the behavior 
of the whole aggregate. The following lemma says that, with high probability, a constant fraction 
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of the aggregate is absorbed in a wide enough (yet still very fine) annulus. Recall the hypothesis 
about the initial cluster in Theorem 9.2. 


Lemma 9.5. [12, Lemma 6] Let I be as above. There exist e,e\,e 2 > 0 and p < 1 such that for all 
n large enough, for all n 1 / 3 < k < en 2 and x\, ...,x k E U n n{U\Ug}, and for all S C [U n n {U\Ug}] 
where (3 £ I. Moreover as in the previous lemma, let B be a set which intersects at most e± fraction 
of the ILJW ’s between Ug and U. k i /2 , then, 

n 


P(|A 


U B; xi ,..., x k -A- U n PI {U\U^ + 



\SUB\ < e 2 k) < p k 


The above lemma says that if at any time the aggregate is a subset of U\Ug along with some e\ 

k i /2 , then in the next k rounds a constant fraction of 

n 

particles fill up sites in U\U k i /2 with high probability. 

^ ' n. 


fraction of the U'^’s between Ug and U 


P+ 


Before proving the above lemmas we prove Theorem 9.2. 

Proof of Theorem 9.2. Recall from (1.6) that Ug = U( Q ). And clearly U( Q ._c v /«) C Ug+c? iv ^ f° r 
some Ci depending on C, (due to the same reasoning as in (4.4)). Thus it suffices to prove that 
the growth cluster will not touch any new points in C Ug+c'i f° r some Ci- 

The proof is along the lines of [ ] and consists of inductively constructing a sequence of aggregates 

A(j) by pausing the particles at different distances mj from the starting set V 2 ,n- If kj is the number 
of paused particles, we choose the next distance rrij+i , at which we pause the particles again, in 
terms of mj and kj. We iterate this procedure until there are less than n 1 /' 3 paused particles. Since 
there are too few particles now, we naively bound the growth of the cluster from this point on. 
Define A^,mj, Pj, kj as follows: A( 0 ) b e th® initial set of filled sites as in the hypothesis of Theorem 
9.2. m\ = \fen and be the cluster after en 2 random walks were released uniformly from IL^n 
(see Setup 1) and stopped on hitting Ug, where fii = /3 + Correspondingly let Pi be the the 
sequence of stopped particles and k\ = |Pi|. We now define , Pj for all j > 1: For j > 1 define 
mj+i in the following way: (this is where we make the adaptations from [ ] since in our setting 

the initial cluster has some particles and hence not all the sites in U n n Ug^ are empty. Consider 
the sequence 

m j + + 2 k 1 ^ 2 , mj + 3 k 1 ^ 2 , - 

Recall U® from (4.5). Let l be the smallest integer such that at most ei’th fraction of the ‘shells’ 
IJi* 3 between U , , 1 / 2 ., and U , 1 / 2 ,, contain particles from the initial cluster A( 0 -) 

p-\-\plj-\-yt J/^ p~r[Ttlj-\-ZKj i/Tl ' ' 

1 /2 

(ei is the same as in the hypothesis of Lemma 9.5). We define and 1 to be mj + (£ — 1) k -' 
1/2 

and mj + £kj 7 respectively. That is we define nij+x to be the smallest number of the right scale 
such that there are enough empty sites between Ug +Tn / n and Ug+ m ',/n. 

Let Aq- +1 ) = A{A^y, Pj —> U n n {U\Ug+ mj - + i/n}); an d Pj+ 1 = Pj n {U\Ug +mj . + i/ n }) 

and let kj +1 = \Pj+i \ ■ Also let J be the (random) first time at which kj < n 1//3 . We do not 
stop the particles after this point and hence by construction, A^p = A(j+i) for any j > .7+1. 
As mentioned earlier the Abelian property guarantees that has the same law as the IDLA 

cluster with m 0 = en 2 particles. Naively A( J+1 ) C U n n {U\Ug +mj / n+ c(e'n+n 1 /3)/ n } for some 
constant C = C(U),. This follows by the following crude bound: since there are at most n 1 ^ 3 
particles left and initially at most e'n many shells U® filled up, the growth from J onwards cannot 
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be more than by e'n + n 1 / 3 shells U®. By Lemma 9.5 for some €2 < 1, 

P[3 2 < j < J A n : kj > (1 — € 2 )^ 1 ] < P[3 2 < j < J A n : kj > (1 — £ 2 )^- 1 ] < np n 1 . 

Since k\ < en 2 by hypothesis this implies that with probability at least 1 — np” 1/3 , we have J < n 
and kj decreases geometrically at each step, hence, 

mj = m\ + (mi - mi) + k\^ + [ml 2 ~ m 2 ) + ... + {m'j_i - mj-i) + k l J\ (9.7) 

1 J_1 

- mi + k ^' \- (t-g^i/a + B m ^ 

Now we finish by bounding Y2e=i( m e ~~ m f)- Notice that by construction at least ei’th fraction 
of all the shells between m } — rrif intersect the initial cluster. Also by hypothesis on the initial 
cluster number of such shells is at most e'n < en/D (recall the hypothesis of Theorem 9.2, we 
specify the value of D here). Thus ei Yle=i( m e ~ m — T) n anc ^ h ence choosing D = we get 
that { m 't ~ m () < en. Also naively ki < en 2 and hence putting everything together with 

probability at least 1 — np n , mj < Cmi = C^/en, for some absolute constant C. Thus ^4(,/+i) 
does not intersect Up +Cl ^ and hence we are done. □ 

We now provide the proofs of Lemmas 9.4 and 9.5. 

Proof of Lemma 9.f We stop the random walk f when it hits U^ +r / 2 n (call the point at which it 
stops as y) and look at the ball of radius cr/n (for some small constant c) around it. Because of the 
smoothness assumption on U even if y lies close to the boundary of U, \{B(y, cr/n) H U n }| > 2 dr 2 , 
for some d = d(c,U). Now if e\ < d then \{B(y,cr/n) fl U n }\5| > dr 2 . Also by hypothesis there 
are at most e\r shells U® initially with particles at time 0 and hence \{B(y, cr/n) n U n }\{S' U 
initial particles}| > dr 2 / 2 by choosing suitably small ei = ei(c).The result now follows from the 
heat kernel estimate in [L6, Theorem 5.2] which shows that the probability that the random walk 
started at y hits the set { B(y,cr/n) n U ri }\{5 U initial particles} before exiting B(y,r/2n) is a 
constant independent of r, n. □ 

Proof of Lemma 9.5 The proof of [12, Lemma 6] uses [12, Lemma 5]. The same arguments as in 
[12, Lemma 6] using Lemma 9.4 above instead of [12, Lemma 5] completes the proof. □ 
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